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Abstract. We review some recent attempts to extract information about 
the natm'e of quantum gravity, with and without matter, by quantum field 
theoretical methods. More specifically, we work within a covariant lattice 
5^ I approach where the individual space-time geometries are constructed from 

fundamental simplicial building blocks, and the path integral over geome- 
tries is approximated by summing over a class of piece-wise linear geome- 
tries. This method of "dynamical triangulations" is very powerful in 2d, 
where the regularized theory can be solved explicitly, and gives us more in- 
sights into the quantum nature of 2d space-time than continuum methods 
are presently able to provide. It also allows us to establish an explicit re- 
lation between the Lorentzian- and Euclidean-signature quantum theories. 
Analogous regularized gravitational models can be set up in higher dimen- 
sions. Some analytic tools exist to study their state sums, but, unlike in 
2d, no complete analytic solutions have yet been constructed. However, a 
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great advantage of our approach is the fact that it is well-suited for numer- 
ical simulations. In the second part of this review we describe the relevant 
Monte Carlo techniques, as well as some of the physical results that have 
been obtained from the simulations of Euclidean gravity. We also explain 
why the Lorentzian version of dynamical triangulations is a promising can- 
didate for a non-perturbative theory of quantum gravity. 



1. Introduction 

There is at present no satisfactory theory of four-dimensional quantum 
gravity. This is partly related to the conceptual questions that arise when 
dealing with fluctuating geometries without reference to any background 
metric. Addressing these seems to call for a genuinely non-perturbative 
formulation of quantum gravity. The need for a non-perturbative approach 
persists in formulations where gravity appears embedded into a larger the- 
ory, such as string or M-theory. Although there have been attempts to 
identify non-perturbative structures within these unified theories (sec [1, 2] 
and many successive papers), they so far seem to have raised as many 
questions as they have been able to answer. In these lectures we will review 
an alternative non-perturbative approach to quantum gravity. It is more 
conservative in spirit, in that it does not conjecture the existence of any 
radically new physical principles or symmetries. This so-called "dynamical 
triangulations" approach^ can be formulated entirely within the framework 
of ordinary quantum field theory, while taking into account that the govern- 
ing symmetry principle is reparametrization- (diffeomorphism-) invariance. 

Our ansatz can be made in any dimension, and has been solved ana- 
lytically in two space-time dimensions. In higher dimensions one so far has 
to rely largely on numerical simulations, and up to now only the case of 
Euclidean signature has been studied in detail. Although gravity in two 
dimensions has no genuine dynamical degrees of freedom, it does provide a 
testing ground for addressing some of the conceptual problems of quantum 
gravity explicitly, for example, how to define a notion of reparametrization- 
invariant length in a fluctuating geometry or how to define correlation func- 
tions without reference to a background geometry. In 2d Euclidean gravity, 

'^The method of dynamical triangulations was introduced in the context of string the- 
ory and 2d quantum gravity in [3, 4, 5], and subsequently extended to higher-dimensional 
Euclidean quantum gravity [6, 7]. An extensive review covering the developments up to 
1996 can be found in the book [8]. A more recent summary is contained in [9], while the 
review [10] deals with a variety of lattice approaches to four-dimensional quantum grav- 
ity, including dynamical triangulations. The use of dynamical-triangulations methods in 
Lorentzian gravity was pioneered in [11, 12, 13]. 
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a number of questions can be answered by continuum functional methods 
(within Liouvihe quantum field theory), but some of the most interesting 
problems involving geodesic distances can only be addressed in the non- 
pcrturbativc setting of dynamical triangulations. The fact that in this latter 
approach quantum gravity is obtained as the scaling limit of a lattice theory 
shows that, contrary to common belief, discrete methods can be used also 
for reparametrization-invariant theories, as long as the discretization takes 
place directly on the space of geometries (that is, on the space of metrics 
modulo diffeomorphisms). 

Another point that has been analyzed in the two-dimensional model is 
the analytic continuation between geometries of Lorentzian and Euclidean 
signature. Canonical quantization attempts of (Lorentzian) gravity usually 
take as their starting point globally hyperbolic manifolds equipped with a 
non-degenerate Lorentzian metric (which gives rise to a causal structure). 
However, it is a priori unclear to what extent these essentially classical 
structures should be preserved in the quantum theory. For example, in a 
covariant path-integral approach it is not obvious that a causality constraint 
should be imposed on each individual space-time configuration contributing 
to the amplitude. This idea was first advocated by Teitelboim [14], and 
more recently in [15, 16]. One argument in its favour is that it seems hard 
to imagine how any notion of causality could emerge in the full quantum 
theory unless it had been imposed in some form on the individual histories 
in the first place. 

Two-dimensional quantum gravity is an ideal testing ground for such 
ideas, since it can be solved explicitly. We will show that by restricting 
the state sum to discrete geometries with a causal structure one obtains a 
theory of quantum gravity with the following features [11, 12]. 

(1) The expectation values of reparametrization-invariant space-time dis- 
tances have canonical dimensions. In other words, in spite of strong 
fluctuations in the geometry, there is still a sense in which the quan- 
tum space-time is two-dimensional. 

(2) When matter with conformal charge c < 1 is coupled to this Lorentz- 
ian quantum gravity theory, both property (1) and the properties of 
the conformal field theory (critical exponents etc.) remain unchanged. 

There is an alternative way of arriving at the same theory, namely, 
by starting from the geometric configurations of the 2d Euclidean gravity 
theory, and removing its "baby universes" in a systematic manner. (The 
proliferation of these branching structures is responsible for most of the 
fractal and geometric properties of the Euclidean theory.) In this way one 
creates a many-to-1 correspondence between Euclidean and Lorentzian ge- 
ometries. Lorentzian 2d quantum gravity appears in this construction as a 
"renormalized" version of the Euclidean quantum theory [17]. 
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Loosening the requirement of causality for the individual space-time his- 
tories, one may allow for changes in the topology of the spatial (constant- 
time) slices as a function of time, without changing the overall topology of 
space-time. This is of course tantamount to reintroducing geometries with 
baby universes. We will show below that once this process is permitted, 
it will totally dominate the structure of space-time. The resulting theory 
(which is equivalent to Euclidean Liouville quantum gravity) has the fol- 
lowing features. 

(1) The expectation values of geodesic space-time distances have anoma- 
lous dimensions, and the intrinsic HausdorfF dimension of quantum 
space-time is four, and not two [18, 19]. 

(2) When matter with conformal charge c < 1 is coupled to this Euclidean 
quantum gravity theory, both the fractal properties of the gravita- 
tional sector and the critical exponents of the conformal field theory 
are changed. 

We therefore may say that the interaction between gravity and matter is 

weak in the Lorentzian 2d gravity model, at least as long as c < 1. Allowing 
for baby-universe creation leads to a strong coupling between matter and 
gravity: the fractal properties of space-time become a function of the matter 
content and in turn the back-reaction of the fluctuating geometry changes 
the critical properties of the matter. The coupling of Ising models to 2d 
gravity provides a particularly clear geometric illustration of the role played 
by the baby universes [20] . 

If one transcends the so-called c = 1 barrier, even Lorentzian gravity ex- 
hibits a strong gravity-matter interaction, leading to a change in the fractal 
properties of space-time [13]. In the Euclidean case (i.e. including baby uni- 
verses), it has been known for a long time that beyond c = 1 the space-time 
disintegrates completely into fractals. These so-called branched polymers 
can be viewed as trees of baby universes of cut-off size. However, this is not 
what happens in the Lorentzian case. One does observe a phase transition 
for c > 1 (to be precise, the transition must take place somewhere in the 
interval 1 < c < 4 [13, 21, 22]), but the new geometry is less pathological 
than the branched polymers. In particular, the critical matter exponents 
retain their Onsager values. 

Many of the results quoted above for c > 1 have been obtained by Monte 
Carlo simulations. Numerical methods have been very successful in the 
study of discretized quantum gravity. Moreover, in two-dimensional quan- 
tum gravity there has been a fruitful interaction between theory and "exper- 
iment" . Non-rigorous calculations have been verified or falsified by numer- 
ical "experiments", and "observations" have inspired theoretical progress. 
In the Lorentzian case, visualizations of the computer-generated space-time 
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geometries have been useful in understanding the influence of matter on ge- 
ometry [23]. 

In higher dimensions, only partial results have been obtained by ana- 
lytical methods. These include asymptotic estimates for the partition func- 
tions of dynamically triangulated Euclidean gravity in 3d and 4d [24, 25], 
and qualitative descriptions of the extreme branched-polymer and crum- 
pled phases [25, 26, 27]. However, as already mentioned, our regularized 
quantum gravity models can readily be studied by means of computer sim- 
ulations. In a first step, one determines the phase diagram of the discrctized 
theory, in order to locate potential critical points where a continuum limit 
can be obtained. One can then use standard finite-size tools to study the 
scaling of various observables. This method has been very successful and 
has enabled us to study both the fractal properties of space-time, and the 
critical exponents of matter coupled to fluctuating geometries. 

The rest of this paper is organized as follows. In Section 2.1 we describe 
and solve the simplest discretized Lorentzian model of 2d gravity, using a 
"Lorentzian" version of the method of dynamical triangulations. In Section 
2.2, the corresponding continuum limit is obtained, while in Section 3 we 
demonstrate how the inclusion of baby universes changes the structure of 
quantum space-time. A non-perturbative definition of Euclidean quantum 
gravity for arbitrary dimension d by means of dynamical triangulations is 
given in Section 4, where we also discuss the inclusion of matter fields. 
Section 5 outlines the principles for numerical simulations of the model, 
with special emphasis on 2d gravity. In Subsection 5.2, wc describe various 
possibilities of defining notions of (fractal) "dimensionality" in the frame- 
work of Euclidean quantum gravity. We discuss why these "observables" 
are particularly well-suited for use in numerical calculations. Subsection 
5.3 provides the interpretation of the numerical results obtained in 2d Eu- 
clidean quantum gravity. Section 6 describes the numerical approach to 
higher-dimensional Euclidean quantum gravity. Finally, in Section 7, we 
outline the future perspectives for both Euclidean and Lorentzian lattice 
gravity. 

2. Lorentzian gravity in 2d 

2.1. THE DISCRETE MODEL 

In order to solve two-dimensional quantum gravity in a path-integral formu- 
lation, one has to "count" geometries. Since for a fixed space-time topology 
the gravitational action consists only of the cosmological term, all geome- 
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tries of a fixed space-time volume contribute with the same weight^. In 
the simplest model of two-dimensional gravity, no spatial topology changes 
are allowed either. For simplicity, we choose the spatial slices (at constant 
"time") to be circles, so that the overall topology of space-time is of the 
form X [0, 1] (for other spatial boundary conditions, see [28]). This mim- 
icks the situation in classical gravity, where one usually works with globally 
hyperbolic space-times [29] . We will use a natural class of triangulations of 
the cylinder x [0,1], to which we will assign edge lengths (i.e. a dis- 
cretized metric) in such a way that each two-dimensional geometry carries 
a discrete causal structure. Our simplicial space-times will have a preferred 
foliation into a discrete set of circular slices, consisting only of space-like 
edges. The foliation parameter t can be interpreted as a discrete version of 
"proper time" . The main task of our dynamical-triangulations approach is 
to count the geometries in this class, that is, to perform the state sum in 
the regularized context, and then attempt to take a continuum limit. 

The geometry of each spatial slice is uniquely characterized by the length 
assigned to it. In the discretized version, the length L will be "quantized" 
in units of a lattice spacing a, i.e. L = I ■ a where / is an integer. A slice 
will thus be defined by / vertices and / links connecting them. To obtain a 
2d geometry, we will evolve this spatial loop in discrete steps. This leads 
to a preferred notion of (discrete) "time" t, where each loop represents a 
slice of constant t. The propagation from time-slice t to time-slice t + 1 is 
governed by the following rule: each vertex i at time t is connected to ki 
vertices at time t + 1, ki > 1, by links which are assigned squared edge 
lengths — a^. The ki vertices, ki > 1, at time-slice t + 1 will be connected 
by fej — 1 consecutive space-like links, thus forming k^ — 1 triangles. Finally 
the right boundary vertex in the set of ki vertices will be identified with 
the left boundary vertex of the set of /cj+i vertices. In this way we get a 
total of J2i=i{ki ~ 1) vertices (and also links) at time-slice t + 1 and the two 
spatial slices are connected by Yl\=i h = k + k+i triangles (see Fig. 1). 

The elementary building blocks of a geometry are therefore triangles 
with one space- and two time-like edges. We define them to be flat in the 
interior. A consistent way of assigning interior angles to such Minkowskian 
triangles is described in [30] . The angle between two time- like edges is jtt = 
— arccos | , and between a space- and a time-like edge 7st = f + ^ arccos | , 
summing up to 7^ + 27^^ = it. The sum over all angles around a vertex with 
j incoming and k outgoing time-like edges (by definition j, > 1) is given 
by 2-17 + {A — j — k) arccos | . The regular triangulation of flat Minkowski 
space corresponds to j = A; = 2 at all vertices. The volume of a single 
triangle is given by ■ 

■^Maybe surprisingly, in the framework of dynamical triangulations also higher- 
dimensional gravity can be reduced to a pure counting problem, c.f. Section 4.3. 




One may view these geometries as a subclass of all possible triangula- 
tions that allow for the introduction of a causal structure. Namely, if we 
think of all time-like links as being future-directed, a vertex v' lies in the 
future of a vertex v iff there is an oriented sequence of time-like links lead- 
ing from V to v' . Two arbitrary vertices may or may not be causally related 
in this way. 

In quantum gravity we are instructed to sum over all geometries con- 
necting, say, two spatial boundaries of length Li and L2, with the weight 
of each geometry g given by 

e^^l^l, S[g\=Kj./^^ (in 2d), (1) 

where A is the cosmological constant. In our discretized model the bound- 
aries will be characterized by integers li and Z2, the number of vertices or 
links at the two boundaries. The path-integral amplitude for the propaga- 
tion from geometry li to I2 will be the sum over all interpolating surfaces 
of the kind described above, with a weight given by the discretized version 
of (1). Let us call the corresponding amplitude G^^\l\,l2)- We thus have 

00 

G^i\hM) = Y.Gi\h,l2;t), (2) 

t=i 

00 

G^^\h,h;t) = J2Gx\h,l;l)lG^^\l,l2,t-l), (3) 
1=1 

G^\h,h;l) = f E e^'"'^^'^ (4) 

^ {ki,...,kii} 

where A denotes the bare cosmological constant^ (we have absorbed the 
finite triangle volume factor), and where t denotes the total number of 
time-slices connecting li and l2- 



^One obtains the renormalized (continuum) cosmological constant A in (1) by an 
additive renormalization, see below. 
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From a combinatorial point of view it is convenient to mark a vertex on 
the entrance loop in order to get rid of the factors I and 1/Z in (3) and (4), 
that is, 

Gx{li,l2;t) = hG'^^\lul2;t) (5) 

(the unmarking of a point may be thought of as the factoring out by (dis- 
crete) spatial diffeomorphisms) . Note that Gx{li,l2',^) plays the role of a 
transfer matrix, satisfying 

Gx{li,l2;h + t2) = ^Gx{li,l;h) Gx{l,l2;t2) (6) 

I 

Gx{h,l2;t + 1) = ^Gx{li,l;l)Gx{l,l2;t). (7) 
I 

Knowing Gxihjh', 1) allows us to find Gx{h,l2',t) by iterating (7) t times. 
This program is conveniently carried out by introducing the generating 
function for the numbers Gxihjh'jt), 

Gx{x,y;t) = Y,^'y'Gx{k,l;t), (8) 

k,l 

which we can use to rewrite (6) as 

/dz 
^—Gx(,x,z-^;h)Gx(,z,y;t2), (9) 

where the contour should be chosen to include the singularities in the com- 
plex z-plane of Gxix, z^^; ti) but not those of Gx{z, y; t2). 

One can either view the introduction of Gx{x, y; t) as a purely technical 
device or take x and y as boundary cosmological constants, 

X = e'^''', y = e*^°", (10) 

such that x'' = e^-^i"-^ becomes a boundary cosmological term, and similarly 
for = 6*'^°'*'. Let us for notational convenience define 

g = e'^-" (11) 

(not to be confused with the symbol for the continuum metric). For the 
technical purpose of counting we view x, y and g as variables in the complex 
plane. In general the function 

G{x,y-g-t) = Gx{x,y-t) (12) 

will be analytic in a neighbourhood of (x, y, g) = (0, 0, 0). 
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From the definitions (4) and (5) it follows by standard techniques of 
generating functions that we may associate a factor g with each triangle, 
a factor x with each vertex on the entrance loop and a factor y with each 
vertex on the exit loop, leading to 

(13) 

Formula (13) is simply a book-keeping device for all possible ways of evolv- 
ing from an entrance loop of any length in one step to an exit loop of any 
length. The subtraction of the term 1/(1 — gx) has been performed to ex- 
clude the degenerate cases where either the entrance or the exit loop is of 
length zero. 

Prom (13) and eq. (9), with ti = 1, we obtain 

G{x,y;g;t) = G{-^,y; g;t - 1). (14) 

I — gx i- — gx 

This equation can be iterated and the solution written as 

(15) 

where Ft{x) is defined iteratively by 

Let F denote the fixed point of this iterative equation. By standard tech- 
niques one readily obtains 



*^ ^ l-xF + F^t+^{x-Fy ^~ 2g ' ^^^^ 
Inserting (17) in eq. (15), we can write 



F2*(1 - F- 



2\2 



Gix, y; g, t) = ; .1 , ^ (18) 

,(19) 



{At - Btx)iAt - Bt{x + y) + Ctxy) 



{l-xF)- F2*+i {F-x) {l-xF){l-yF)- F"^^ {F-x){F-y) 

where the time-dependent coefficients are given by 

y4t = l-F2*+2, Bt = F{1- F^^), Ct = F2(1 - f2*-2). (20) 
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The combined region of convergence to the expansion in powers g^x^y"^, 
vahd for all t is 

\9\<\. W<1, H<1- (21) 



2.2. THE CONTINUUM LIMIT 

The path integral formalism we are using here is very similar to the one used 
to represent the free particle as a sum over paths. Also there one performs a 
summation over geometric objects (the paths), and the path integral itself 
serves as the propagator. From the particle case it is known that the bare 
mass undergoes an additive renormalization (even for the free particle), 
and that the bare propagator is subject to a wave-function renormalization 
(see [8] for a review). The same is true in two-dimensional Euclidean grav- 
ity, treated in the formalism of dynamical triangulations [8]. The coupling 
constants with positive mass dimension, i.e. the cosmological constant and 
the boundary cosmological constants, undergo an additive renormalization, 
while the partition function itself (i.e. the Hartle-Hawking-like wave func- 
tion) undergoes a multiplicative wave-function renormalization. We there- 
fore expect the bare coupling constants A, Aj and Aq to behave as 

A = %+A, Xi = 9^ + X, A„ = ^ + r, (22) 

where k^X,Y denote the renormalized cosmological and boundary cosmo- 
logical constants. If we introduce the notation 

gc = e'^\ Xe = e^^^S yc = e''^'% (23) 

for critical values of the coupling constants, it follows from (10) and (11) 
that 

5 = 5ce^"'^, x = xce-^, y = 2/ce^"^. (24) 

The wave-function renormalization will appear as a multiplicative cut-off 
dependent factor in front of the bare "Green's function" G{x,y;g;t), 

G^{X,Y;T) = lim a" Gix,y -git), (25) 

a— >0 

where T = at, and where the critical exponent t] should be chosen so that 
the right-hand side of eq. (25) exists. In general this will only be possible 
for particular choices of gc,Xc and yc in (25). 

The basic relation (6) can survive the limit (25) only if ry = 1, since 
we have assumed that the boundary lengths Li and L2 have canonical 
dimensions and satisfy = ali. 
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A closer analysis reveals that only at gc = 1 /2 one can obtain a sensible 
continuum limit. It corresponds to a purely imaginary bare cosmological 
constant A^. := Cxjc? = —iXMljc?. If we want to approach this point from 
the region in the complex g'-plane where G(x, y; g\ t) converges it is natural 
to choose the renormalized coupling A imaginary as well, A = iA, i.e. 

Ini 

X = i ^ + iA. (26) 
a"' 

One obtains a well-defined scaling limit (corresponding to A G IR) by letting 
A — > Ac along the imaginary axis. The Lorentzian form for the continuum 
propagator is obtained by an analytic continuation A — >■ —iA in the renor- 
malized coupling of the resulting Euclidean expressions. 

At this stage it may seem that we are surreptitiously reverting to a 
fully Euclidean model. We could of course cquivalcntly have conducted the 
entire discussion up to this point in the "Euclidean sector", by omitting 
the factor of —i in the exponential (1) of the action, choosing A positive 
real and taking all edge lengths equal to 1. However, from a purely Eu- 
clidean point of view there would not have been any reason for restricting 
the state sum to a subclass of geometries admitting a causal structure. The 
associated preferred notion of a discrete "time" allows us to define an "an- 
alytic continuation in time". Because of the simple form of the action in 
two dimensions, the rotation 

J dx dt\/— del gior i J dx dteuVdetg^ (27) 

to Euclidean metrics in our model is equivalent to the analytic continuation 
of the cosmological constant A. 

From (18) or (19) it follows that we can only get macroscopic loops 
in the limit a ^ if we simultaneously take a;,y — > 1. (For gc = —1/2, 
one needs to take x,y — *■ —1. The continuum expressions one obtains are 
identical to those for gc = 1/2.) Again the critical points correspond to 
purely imaginary bare boundary cosmological coupling constants. We will 
allow for such imaginary couplings and thus approach the critical point 
\ = = ^ from the region of convergence of G{x,y; g;t), i.e. via real, 
positive X, Y where 

Xi = iXa, Ao = iYa. (28) 

Again X and Y have an obvious interpretation as positive boundary cosmo- 
logical constants in a Euclidean theory, which may be analytically continued 
to imaginary values to reach the Lorentzian sector. 
Summarizing, we have 

5 = ^e-^"'^^(l-Aa'), (i.e. F = l-aVA) (29) 
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as well as 

X = e-^" ^l-aX, y = e""^ ^ 1 - aY, (30) 

where the arrows — in (29) and (30) should be viewed as analytic coupling 
constant redefinitions of A, X and Y, which we have performed to get rid 
of factors of 1/2 etc. in the formulas below. With the definitions (29) and 
(30) it is straightforward to perform the continuum limit of G{x,y;g;t) as 
ix,y,9) {xc,yc,9c) = (1,1,1/2), yielding 

4A e-2^^ 

Ga(X,Y;T) = —— — j= 

^ ^ (VA + X)+e-2VAT(^_X) 

^ (31) 



( + x)( Va + y) - e-2VAr(y^ - X) ( - y) ■ 

For T ^ oo one finds 

GAiX, Y; T) ^ ^, (32) 

^ ' ' ^ (X + VA)2(y + VA) ^ ' 

while the limit for T — ^ is 

GA(X,y;r = 0) = ^^, (33) 

in accordance with the expectation that 

GA(Li,L2;r = 0) = <^(Li-L2). (34) 

The general expression for G\{Li, L2; T) can be computed as the inverse 
Laplace transform of formula (31), yielding 

''^^''^'''^'^^ = — S^wx? i j ' ^''^ 

where Ii{x) is a modified Bessel function of the first kind. The asymmetry 
between Li and L2 arises because the entrance loop has a marked point, 
whereas the exit loop has not. The amplitude with both loops marked 
is obtained by multiplying with L2, while the amplitude with no marked 
loops is obtained after dividing (35) by Li. Quite remarkably, our highly 
non-trivial expression (35) agrees with the loop propagator obtained from 
a bona-fide continuum calculation in proper-time gauge of pure 2d gravity 
by Nakayama [31]. 

The basic result (31) for Ga{X,Y]T) can be derived by taking the 
continuum limit of the recursion relation (14). By inserting (29) and (30) 
in eq. (14) and expanding to first order in the lattice spacing a we obtain 

Ag^(X, Y-T) + ^ [{X^ - A)Ga(X, y ; T)] = 0. (36) 
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This is a standard first-order partial differential equation which should be 
solved with the boundary condition (33) at T = 0, since this expresses the 
natural condition (34) on Ga(Li,L2). The solution is thus 

c^(x,r;r, = ^!§^j^^^, (37) 

where X{T;X) is the solution to the characteristic equation 

^ = _(X2-A), X{T = 0)=X. (38) 
It is readily seen that the solution is indeed given by (31) since we obtain 

X(T) = VA'f (39) 
(VA + X)+e-2VAT(^_X) ^ ^ 

If we interpret the propagator Ga(Li, L2; T) as the matrix element be- 
tween two boundary states of a Hamiltonian evolution in "time" T, 

GA(Li,L2;r) = (Li|e-^^|L2) (40) 

we can, after an inverse Laplace transformation, read off the functional 
form of the Hamiltonian operator H from (36), 

The corresponding Hamiltonian for the propagator of unmarked loops 
is given by 

The Hamiltonian (42) is Hermitian with respect to the natural mea- 
sure LdL, which has its origin in the basic completeness relation (3) for 
the transfer matrix with unmarked entrance and exit loops. If one wants 
to construct a unitary evolution with respect to the "time" -parameter T 
appearing in the transfer-matrix approach, one can simply exponentiate 
{iHT). 

However, we should point out that there is an alternative to the analytic 
continuation T — > —iT if one wants to relate the Euclidean and Lorentzian 
sectors of the theory. The combination s/AT appearing as an argument in 
(35) arises in taking the continuum limit of powers of the form in expres- 
sions like (18), (19), where F is defined in (17). There are two aspects to a 
possible analytic continuation of F*. The power t in should clearly not be 
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continued, since it is simply an integer counting the number of iterations of 
the transfer matrix. However, the function F itself does refer to the action, 
because the dimcnsionlcss coupling constant g = c*'^"*"' is the action for a 
single Lorentzian triangle. (For added clarity we have distinguished between 
the lattice spacings in time- and space-directions, and called them at and 
a;.) From the expression for F in terms of g in (17), wc have F = 1 — ^/0^0; A. 
The analytic continuation of F in time, from Euclidean to Lorentzian time, 
corresponds to the substitution at — > —i at under the square-root sign, and 
thus becomes equivalent to the continuation A — ^ —iK in the cosmologi- 
cal constant, as already remarked below cq. (27). The subtleties associated 
with the analytic continuation in the "time" -parameter T appearing in 
a transfer-matrix formulation of quantum gravity were first discussed in 
[32, 33] in the context of a square-root action formulation. Similar difficul- 
ties will also be present in higher-dimensional gravity, where the analytic 
continuation from Euclidean metrics to Lorentzian metrics cannot be ab- 
sorbed by a continuation in A alone. To conclude, it is not obvious how to 
choose the correct analytic continuation back to Lorentzian signature, once 
the continuum limit has been taken. The continuation T ~^ —iT leads to 
a unitary theory. For the continuation A — > —iA, we have not yet found a 
scalar product which makes the corresponding evolution operator unitary. 

3. Topology changes and Euclidean quantum gravity 
3.1. BABY UNIVERSE CREATION 

In our non-perturbative regularization of 2d quantum gravity wc have so 
far not included the possibility of topology changes of space. Wc will now 
show that if one allows for spatial topology changes, one is led in an essen- 
tially unambiguous manner to a different theory of two-dimensional quan- 
tum gravity, where space-time is much more fractal, and which agrees with 
Euclidean 2d quantum gravity as defined by Euclidean dynamical triangu- 
lations or Liouville theory. 

By a topology change of space in our Lorentzian setting we have in mind 
the following: a baby universe may branch off at some time T and develop 
in the future, where it will eventually disappear in the vacuum, but it is 
not allowed to rejoin the "parent" universe and thus change the overall 
topology of the two-dimensional manifold. This is a restriction we impose 
to be able to compare with the analogous calculation in usual 2d Euclidean 
quantum gravity. 

It is well known that such a branching leads to additional complica- 
tions, compared with the Euclidean situation, in the sense that, in general, 
no continuum Lorentzian metrics which are smooth and non-degenerate ev- 
erywhere can be defined on such space-times (see, for example, [34] and ref- 
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erences therein). These considerations do not affect the cosmological term in 
the action, but lead potentially to contributions from the Einstein-Hilbert 
term at the singular points where a branching or pinching occurs. 

We have so far ignored the curvature term in the action since it gives 
merely a constant contribution in the absence of topology change. We will 
continue to do so in the slightly generalized setting just introduced. The 
continuum results of [34] suggest that the contributions from the two sin- 
gular points associated with each branching of a baby universe (one at the 
branching point and one at the tip of the baby universe where it contracts 
to a point) cancel in the action. The physical geometry of these configu- 
rations may seem slightly contrived, but they may well be important in 
the quantum theory of gravity and deserve further study. However, for the 
moment our main motivation for introducing them is to make contact with 
the usual non-pcrturbative Euclidean path-integral results. 

We will use the rest of this section to demonstrate the following: once 
we allow for spatial topology changes, 

(1) this process completely dominates and changes the critical behaviour 
of the discretized theory, and 

(2) the disc amplitude W\{L) (the Hartle-Hawking wave function) is uni- 
quely determined, almost without any calculations. 

Our starting point will be the discretized model introduced in Sec. 2.1. In 

this model wc do not directly have a disc amplitude like the Hartle-Hawking 
wave functional. However, as discussed at the beginning of this section, the 
degeneracy of the metric at a point in the interior of the disc is always 
compensated (in the sense of complex contributions to the action) by the 
degeneracy of the metric at the point where the baby universe branched 
off. We can thus define the disc amplitude in Lorentzian gravity as 

w'-'^\x,g) ■.= Y,G^'\x,l2 = l;9;t)=G'-''\x,l2 = l;9), (43) 
t 

and the continuum version as 

{X) =JdT {X, L2 = 0; T) = {X, = 0) (44) 

where the superscript ^''^ indicates the "bare" Lorentzian model without 
spatial topology changes. It follows that 

w^^\x,g)^a-^w!l'\x). (45) 

There are a number of ways to implement the creation of baby universes, 
some more natural than others, but they all agree in the continuum limit, 
as will be clear from the general arguments provided below. Here we discuss 
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Figure 2. A "baby universe" created by a global pinching. 



only the simplest way of implementing such a change. This is shown in Fig. 
2: stepping forward from t to t + 1 from a loop of length li we create a 
baby universe of length Z < Zi by pinching it off non-locally from the main 
branch. 

Accounting for the new possibilities of evolution in each step according 
to fig. 2, the new and old transfer matrices are related by 

h-i 

GxihM; 1) = of {h , ^2; 1) + E hw{h-l, g) of {I, h;!). (46) 

1=1 

The factor li in the sum comes from the fact that the "pinching" shown in 
fig. 2 can take place at any of the li vertices. As before, the new transfer 
matrix leads to new amplitudes G\{li,l2',t), satisfying 

Gxih, WM + i2) = E ^a(^i, ^; h)Gx{l, h; ^2), (47) 
I 

and in particular 

Gx{li,l2;t) = J2Gxili,l;l) Gx{l,l2;t-l). (48) 

Performing a (discrete) Laplace transformation of eq. (48) leads to 
G{x,y;g;t) = 



dz 
2m z 



(x, ; 1) {wix- g)G^;^ {x,z~^-l)) 



Giz,y;g;t-1), 



(49) 



or, using the explicit form of the transfer matrix G^^ {x, z; 1), formula (13), 



G{x,y;g;t) 



dw{x,g) d 
— dx ^^'^^^'3)q^ 



gx 



1- gx 



G(-^,y;g;t-l). 

(50) 
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At this point neither the disc amphtude w{x, g) nor G{x, y; g; t) are known. 
We will now show that they are uniquely determined if we assume that the 
boundary length scales canonically with the lattice spacing, L = al, imply- 
ing a renormalized boundary cosmological constant X with the dimension 
of mass, X = Xc{l — aX). In addition we assume that the dimension of the 
renormalized cosmological constant A is canonical too, g = gc{l — ^Aa^). 
Somewhat related arguments have been presented in different settings in 
[35, 36]. 

It follows from relation (47) that we need 

Gxili,l2,t)''-^aGAiLi,L2;T). (51) 

It is important for the following discussion that G\{li, h'jt) cannot contain 
a non-scaling part since from first principles (sub-additivity) it has to decay 
exponentially in t. By a Laplace transformation, using x = Xc{l — aX) in 
the scaling limit, we thus conclude that 

Gx{x,l2,t)'^GA{X,L2,T), (52) 

and further, by a Laplace transformation in L2, 

Gxix, y; t) ^ a-'GAiX, Y; T). (53) 

We will now show that the scaling of w{x,g) is quite restricted too. The 
starting point is a combinatorial identity which the disc amplitude has to 
satisfy. The arguments are valid both for the disc amplitude in Euclidean 
quantum gravity and the disc amplitude we have introduced for our model 
in (43). We will assume the general scaling 

ui{x, g) = Wns{x-, g) + aPW^iX) + less singular terms (54) 

for the disc amplitude. In the case 77 < the first term is considered absent 
(or irrelevant). However, if 77 > a term like Wns will generically be present, 
since any slight redefinition of the coupling constants of the model will 
produce such a term if it was not there from the beginning. 

We will introduce an explicit mark in the bulk of w{x,g) by differenti- 
ating with respect to g. This leads to the combinatorial identity 



9 = E E G{x, I; g; t) I w{l, g), (55) 

t I 

or, after the usual Laplace transform. 



dw{x,g) y-^ f dz _i dw{z,g) 
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T 




Figure 3. Marking a vertex in the bulk of Wa{X). The mark has a distance T from the 
boundary loop, which itself has one marked vertex. 

The situation is illustrated in Fig. 3. A given mark has a distance t {T 
in the continuum) to the entrance loop. In the figure we have drawn all 
points which have the same distance to the entrance loop and which form 
a connected loop. In the bare model these are all the points at distance t. 
In the case where baby universes are allowed (which we have not included 
in the figTirc), there can be many disconnected loops at the same distance. 
Let us assume a general scaling 



for the time variable T in the continuum limit. Above we saw that the 
bare model without baby universe creation corresponded to £ = 1. With 
the generalization (57) we account for the fact that by allowing for baby 
universes we have introduced an explicit asymmetry between the time- and 
space-directions . 

Inserting (54) and (57) into cq. (56) we obtain 



T = aH^ 



£ > 0, 



(57) 



dw. 



'ns 



-2aP 



2 dWK{X) 

dk 




Z-T) 



dw. 



'ns 



- a- 



(58) 



dz 



2c dZ 



where (x,(/) = {xcQc) in the non-singular part. 

From eq. (58) and the requirement e > it follows that the only con- 
sistent choices for r] are 
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1: 77 < 0, i.e. 

dK 2aP 



liTj4Za,(X.-Z;T)l'-^, (59) 



in which case we get £ = 1; and 
2: 1 < 97 < 2. Here formula (58) sphts into the two equations 



" M 2a= dz 



j dT J dZ Ga{X,-Z;T), (60) 



and 



dw 



ns 



dg 



9=9c 



'^[,t[,ZG4X.-Z;T)^?^. ,61) 

J J Zc aZ 



Wc arc led to the conclusion that e = 1/2 and rj = 3/2, which are pre- 
cisely the values found in Euclidean 2d gravity. Let us further remark 
that eq. (60) in this case becomes 

- = const. Ga{X, L2 = 0), (62) 

which differs from (44) by a derivative with respect to the cosmological 
constant. We will explain the reason for this difference below. Finally, 
eq. (61) becomes 

JdT J dZ Ga{X, -Z; T) ^^q^^ = const., (63) 

which will be satisfied automatically if = 3/2 and e = 1/2, as we will 
show below. 

We will now analyze a possible scaling limit of (50), assuming the canon- 
ical scaling x = Xc{l — aX) and g = gd^— ^Aa^). In order that the equation 
have a scaling limit at all, Xc, gc and Wns{xc, Qc) must satisfy two relations 
which can be determined straightforwardly from (50). The remaining con- 
tinuum equation reads 

a'—GA{X,Y;T) = -a — [(X^ - A)Ga(X, Y; T) 



" dX 



Wa{X)Ga{X,Y;T) . (64) 



The first term on the right-hand side of eq. (64) is precisely the one we have 
already encountered in our original model, while the second term is due to 
the creation of baby universes. Clearly the case ry < (in fact r/ < 1) is 
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inconsistent with the presence of the second term, i.e. the creation of baby 
universes. However, since r) < 2, the last term on the right-hand side of (64) 
will always dominate over the first term. Once we allow for the creation of 
baby universes, this process will completely dominate the continuum limit. 
In addition we get e = rj — 1, in agreement with (61). It follows that r] > 1 
and we conclude that e = 1/2, rj = 3/2 are the only possible scaling ex- 
ponents if we allow for the creation of baby universes. These are precisely 
the scaling exponents obtained from two-dimensional Euclidean gravity in 
terms of dynamical triangulations, as we have already remarked. The topol- 
ogy changes of space have induced an anomalous dimension for T. If the 
second term on the right-hand side of (64) had been absent, this would 
have led to e = 1, and the time T scaling in the same way as the spatial 
length L. 

In summary, in the case {r],£) = (3/2, 1/2) eq. (64) leads to the contin- 
uum equation 

— GAiX,Y;T) = -—[WAiX)GA{X,Y-T)], (65) 

which, combined with eq. (62), determines the continuum disc amplitude 
Wa{X). Integrating (65) with respect to T and using that Ga{Li,L2;T = 
0) = (5(Li-L2),i.e. 

Ga{X,L2 = 0;T = 0) = 1, (66) 

we obtain 

= axK^^^M^'^^^^J- ^^^^ 

Since Wa{X) has length dimension -3/2, i.e. Wl{X) = X^F{VA/X), the 
general solution must be of the form 



Wa{X) = yJ-2KX + 62^3 + c^A?l^. (68) 

From the very origin of Wa{X) as the Laplace transform of a disc amplitude 
Wa{L) which is bounded, it follows that Wa{X) has no singularities or cuts 
for Re X > 0. This requirement fixes the constants 6, c in (68) such that 



Ir.(x) = 6(x-^VA)^X + MVA, (69) 

where the constant 5 is determined by the model-dependent constant in 
(62). This expression for the disc amplitude agrees after a rescaling of the 
cosmological constant with the amplitude W^^\x) from 2d Euclidean 
quantum gravity. 



WA{X) = {X-^^^X) \[)cWi. (70) 
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With Wa{X) substituted into (65), the resulting equation is familiar from 
the usual theory of 2d Euclidean quantum gravity, where it has been derived 
in various ways [18, 35, 36], with T playing the role of geodesic distance 
between the initial and final loop. 

Before showing that the anomalous scaling of the proper time - once 
baby universes are allowed - leads to an intrinsic fractal space-time di- 
mension of four (rather than two), let us comment on the difference be- 
tween the equations for the amplitudes (59)-(61) for {r],e) = (—1,1) and 
(77,6) = (3/2, 1/2) respectively. In the first case there are no baby universes 
and eq. (59) entails that only macroscopic loops at a distance T from the 
entrance loop are important (as illustrated in Fig. 3). On the other hand, 
the term dWA{Z)/dZ which describes the presence of these macroscopic 
loops is absent in eq. (60). This is consistent with eq. (62), which shows 
explicitly that the length of the upper loop in Fig. 3 remains at the cut- 
off scale, and therefore can never become macroscopic. It is also consistent 
with the abundance of baby universes, since at any point in space-time the 
probability for creating a little "tip" of cut-off size will dominate. At the 
same time, the right-hand side of cq. (59), that is. cq. (63), will play no 
role when 1 < rj < 2, being simply equal to a constant. This latter property 
is satisfied automatically, as can be seen by using an equation analogous 
to (65) for the exit instead of the entrance loop. Thus eq. (63) becomes 
proportional to 

j^'drA G(X,L2=0;r) = const., (71) 
proving our previous assertion. 



3.2. THE FRACTAL DIMENSION OF EUCLIDEAN 2D GRAVITY 

If we allow for baby universe creation, the fractal structure of space-time is 

determined by (65) and (70), where T is the time separating the entrance 
and exit loops. As already mentioned, these are exactly the equations gov- 
erning Euclidean quantum gravity, if one replaces T by the geodesic distance 
between the two loops. 

One can solve eq. (65) in the same way as its Lorentzian analogue (36) 
was solved by (37)-(39). We have 

Wa{X{T;X)) 1 
^^^^'^'^^ = W^iX) X{T;X)+Y' ^^^^ 

where X{T;X) is the solution to the characteristic equation 

^ = -Wa{X), X{T = 0)=X. (73) 
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This equation can be solved in terms of elementary functions. In particular, 
one finds 

X(r; X = oo) oc VA coth^ (aA^'^T) , a = \\JV^- (74) 

We may now define a two-point function by contracting the entrance loop 
to a point and closing the exit loop by the disc amplitude. This is shown 
in Fig. 3, except that now the entrance loop must be contracted'^. The 
two-point function Ga{T) in Euclidean gravity can be interpreted as the 
(unnormalized) average over all geometries with two marked points which 
are separated by a geodesic distance T. At the discretized level it is defined 
by an equation analogous to (55), but without summing over t, 

G{t;g) = Y,G{h,l;9;t)lw{l,g). (75) 
I 

As in eq. (58), the continuum limit will be dominated by the non-scaling 
part of w{l;g), i.e. by small I, and the continuum two-point function is 
simply given by 

Ga(T) ~Ga(Li=0,L2 = 0;T). (76) 
Using (72) and the Laplace transforms of Ga and Wa we obtain 

^^^^^ = / 2^ / 2^ Wa{X) X{T-X)+Y - ^''^ 
With the help of the characteristic equation (73) this leads to 

d r dx x{T-x) 

^^^^^^dTf2^im^- ^^^^ 

The contour of integration can be deformed to infinity and we obtain from 
(74) that 

This two-point function may be viewed as the partition function for uni- 
verses with two marked points separated by a geodesic distance T. Since 

''One cannot a priori contract both loops, since the two points in the two-point func- 
tion are separated by a geodesic distance T, and there may be many points at distance 
T from the entrance point, as shown in Fig. 3. They will in general form several con- 
nected loops. However, after solving the model, it turns out that one can just contract 
the exit loop and obtain the two-point function. The reason is that w{x,g) contains a 
non- universal part, which again implies that a typical disc amplitude in Fig. 3 will be of 
cut-off size. This is precisely the contents of eq. (60), whose non-scaling part is given by 

dWris/dz\;,=:c^. 
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we wanted to solve Euclidean gravity, this is our final result. Given (79), 
we can calculate the average space-time volume of such a universe, 

(V)^ = T^f[k^'^T^, F(0) > 0. (80) 

The function F is again expressible in terms of elementary functions and 
one finds 

{V)t^T^ for r<l/A^/^ (81) 

This formula shows that the intrinsic fractal dimension of 2d Euclidean 
quantum space-time is four, whereas a similar derivation in the case of 
Lorentzian gravity yields [11] 

(y)y«T2 for T<1/A1/^ (82) 

The discrepancy in dimension can be explained purely in terms of the baby 
universe structure: since at each point of the two-dimensional Lorentzian 
surface a baby universe can branch off, the resulting fractal Euclidean space- 
time has twice the intrinsic dimension. 

4. Euclidean quantum gravity 

4.1. SOME GENERALITIES 

In the last section we arrived at the 2d Euclidean gravity theory through 
an "extension" of the Lorentzian model. Euclidean quantum gravity can of 
course be defined independently, as the quantization of classical gravity on 
the space of all Riemannian metrics (of positive definite signature) instead 
of the space of (indefinite-signature) Lorentzian metrics. In two dimensions, 
Euclidean gravity has a well-defined continuum path-integral formulation. 
Choosing a conformal gauge-fixing leads to the so-called Liouville gravity. 
Certain aspects of this theory can be solved by a bootstrap approach. In 
higher dimensions the path-integral approach to Euclidean quantum gravity 
is problematic, since the Einstein-Hilbert action is unbounded from below. 
There are various ways of tackling this problem: analytically continuing the 
unstable modes, adding stabilizing higher-derivative terms to the action, or 
defining the theory non-perturbatively via a lattice regularization, such that 
the action is bounded for any finite lattice volume. 

One example of the latter is the dynamical-triangulations approach, 
which has the added bonus of being exactly soluble by combinatorial meth- 
ods in two dimensions. Its continuum limit agrees with continuum Liouville 
quantum gravity, wherever the two fornuilations can be compared. We are 
in fact in the unusual situation that the lattice approach can address and 
answer more questions than the continuum methods. We will not review 
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the combinatorial approach here since the results (for pure gravity) were al- 
ready obtained in Section 3, starting from Lorentzian gravity. For a detailed 
description we refer to [8], chapter 4. The generalization of this Euclidean 
lattice path integral to higher dimensions is straightforward and shares two 
virtues with the 2d case: calculating the partition function for gravity is 
again turned into a combinatorial problem, and the model is well-suited 
for numerical simulations. We have by now a good qualitative understand- 
ing of the phase structure of Euclidean dynamically-triangulated gravity 
in d = 3, 4, although complete analytical solutions of the discretized mod- 
els are still missing. However, the combinatorial nature of the partition 
function gives us some hope that progress can be made also in these cases. 

4.2. DYNAMICAL TRIANGULATIONS 

This lattice approach shares many elements with lattice regularizations of 
ordinary quantum field theory. The main difference lies in the fact that 
the geometric degrees of freedom become dynamical and the lattices are 
therefore no longer part of the inert background structure. The geometric 
quantum fluctuations must be taken properly into account when building 
discretized models of matter interacting with quantum gravity. The field- 
theoretical, non-perturbative Euclidean path integral of such a theory takes 
the general form 



where the integration is over equivalence classes of metrics [g^u] and matter 

fields (pi. The action contains a matter part S,natter, depending on a set of 
matter couplings {Pi}, and a purely geometric part, given by the Einstein- 
Hilbert action with a cosmological term. 



In (83) we have omitted possible boundary terms. Except in d = 2, expres- 
sions of the kind (83) have remained formal, due to the absence of a suitable 
diffeomorphism-invariant integration measure. The lattice formulation is an 
attempt to remedy this situation, by using an intermediate regularization 
of the non-perturbative path integral (83) (sec [37, 8] for reviews). Defining 
a discrete regularization consists of several steps: 

• a discretization of the individual metric manifolds, together with a def- 
inition of discretized geometric "observables" , such as lengths, volumes 
and (scalar) curvature. These are necessary for obtaining a discretized 
version of the action, and for analyzing the physical properties of the 
theory in terms of scaling relations and correlation functions. 




(83) 




(84) 
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• a suitable choice of an integration measure on the space of discretized 
geometries (i.e. equivalence classes of metrics), such that the discrete 
path integral converges. 

• a discretization of the matter sector, which will be closely related to 
standard lattice-regularizations in field theory. 

Let us now describe the dynamical-triangulations regularization of Eu- 
clidean quantum gravity. It consists in replacing the d-dimensional Rie- 
mannian metric continuum manifold by a simplicial manifold constructed 
from equilateral d-dimensional simplices of (geodesic) edge length a. (A 
simplex is a point in d = 0, an edge in d = 1, a triangle in d = 2, a tetra- 
hedron in d = 3, etc.) Using Regge's prescription [38], all quantities can 
be expressed as functions of the squared edge lengths. For example, the 
curvature depends on local deficit angles, which in turn are expressible in 
terms of edge lengths. A simplicial complex is obtained by gluing together 
d-simplices pairwise along (d — l)-dimensional faces (which are themselves 
{d— l)-simplices). Since, unlike in Regge calculus, our edge lengths are not 
variable, all d-simplices have the same size, and the total volume of the sim- 
plicial complex is simply proportional to the number of such cells. Each 
d-simplex is built from simplices of lower dimensionality. It contains d+1 0- 
simplices (vertices), {d + l)d/2 1-simplices (hnks) etc. A lower-dimensional 
subsimplcx is in general shared by a number of d-simplices, called the or- 
der of the subsimplcx. A simplicial complex is a simplicial manifold if the 
neighbourhood of any p-simplex (p < d) has the topology of a (d — p — 1)- 
dimensional sphere. Physically the manifold requirement may be viewed as 
a regularity condition at the cut-off scale, which will be convenient to use 
in our construction. The numbers A^^ of (sub-)simplices of dimension k < d 
are not independent, but (due to the regularity requirement) must satisfy 
a set of so-called Dehn-Sommerville relations, namely, 

^^ = i:(-ir'^( -^J (85) 
together with the Euler constraint 

j2i-l)'N, = x. (86) 

fc=0 

For fixed Euler number x and d = 2, all can be expressed as functions 
of the single variable N2, say. For d = 3, 4, two of the Nk are independent. 

A triangulation T together with an assignment of geodesic edge lengths 
and flat simplex interiors may be viewed as a piecewise linear manifold, 
and provides an explicit coordinate-independent representation of a metric 
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manifold. In this way each triangulation corresponds to a unique equiva- 
lence class of metrics (albeit of piecewise-linear, and not of differentiable 
type). 

When it comes to numerical simulations, it is often convenient to assign 
a label to each vertex. From the list of vertex labels for all d-simplices the 
whole manifold can be reconstructed. Since the labels themselves have no 
physical meaning, the labelling introduces a redundancy. Invariancc under 
permutations of the labels may loosely be regarded as a discrete analogue of 
the diffeomorphism invariance of a differentiable manifold. Also the lower- 
dimensional subsimplices are characterized by their vertex labels. More- 
over, the regularity requirement implies that we cannot have two different 
(sub-)simplices with the same set of vertex labels. For each triangulation 
T with Nq vertices the number of different labellings equals No\/C{T), 
where C(T) is the order of the automorphism group of T. We can there- 
fore distinguish between labelled triangulations T and abstract unlabelled 
triangulations T. As mentioned above, different T's (with fixed topology) 
correspond to different equivalence classes of piecewise linear metrics as- 
signed to the underlying manifold, and allow us to work directly with a 
reparameterization-invariant set of geometries. 

Since the simplices are flat on the inside, curvature is located (distri- 
butionally) at simplices of lower dimension. Circulating around a (d — 2)- 
dimensional simplex, the contributing angles will in general not add up to 
27r. The resulting deficit angle depends on the number of simplices meeting 
at the subsimplex. We conclude that the curvature is concentrated at the 
{d — 2)-dimensional subsimplices of the triangulation. The Einstein-Hilbert 
action for a dynamically triangulated manifold in dimension d assumes the 
simple form 

Seh = HdNd - Kd-2Nd-2, (87) 

with the two dimensionless coupling constants and Kd-2- As usual in 
lattice regularizations, the lattice spacing a has disappeared from the for- 
mulation and will have to be reintroduced in the scaling limit. We are using 
the Einstein-Hilbert action because of its simplicity; one could in principle 
consider also the inclusion of higher-order curvature terms. 

With each simplicial lattice described above one can associate a dual 
lattice, whose vertices are located at the centres of the simplices of the 
original lattice. In a similar way we can associate to each p-simplex a dual 
object of dimension d—p. For example, the dual links connect the centres of 
neighbouring simplices, and the dual of a (d — 2)-simplex is a closed loop (a 
two-dimensional object) whose length is equal to the number of d-simplices 
of the original lattice which share the {d — 2)-simplex. Since in d dimensions 
each simplex has d+1 neighbours, the dual lattices have the form of graphs 
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of a scalar + -theory (that is, aU their vertices are (d + l)-valent), but 
with a local d-dimensional topological structure. 

The simplicial structures described above possess a natural notion of 
length for any path connecting two points, since the equivalence class of 
metrics is uniquely determined. To simplify matters, we will only consider 
certain sets of discretized paths on the simplicial manifold. The first set is 
given by paths which run along the links of the original simplicial lattice, 
and the other by paths running along the links of the dual lattice. In either 
case we may define a distance between points on the lattice or its dual as 
the number of edges of the shortest path connecting the two. At first glance 
these definitions seem different from the standard notion of a geodesic dis- 
tance, but they all coincide in the scaling limit, up to trivial numerical 
factors^ . Numerical tests of this assumption will be discussed below. 



4.3. THE FUNCTIONAL INTEGRAL 

The association of triangulations with equivalence classes of metrics mo- 
tivates the use of the discrete sum over d-dimensional triangulations as a 
discretized analogue of the diffeomorphism-invariant integration measure 
in (83), 

^M^E^' (88) 



/ 



where the sum is taken over unlabelled simplicial manifolds. The need for 
including the symmetry factor C(T) becomes apparent when one rewrites 
the right-hand side of (88) as a sum over labelled triangulations T, 

1 1 
^C(T)^^]^^- ^^^^ 

In order that a discretized path integral with this choice of measure lead to 
a theory with a well-defined thermodynamic limit, the number of triangu- 
lations with fixed volume should grow at most exponentially with as 
Nd oo. This is not the case unless we fix the space-time topology (usually 
to that of a sphere S*^); otherwise the growth is factorial. This property has 
been proven for d = 2 and arbitrary topology [39] and for simply connected 
manifolds in d = 3, 4 [24, 25]. 

^The situation is the same as for a regular 2d quadratic lattice in flat space: if we 
are only allowed to connect vertices along the lattice links, the lattice distance between 
different lattice points can differ by as much as a factor \/2 from the Euclidean distance 
in flat space. However, in the scaling limit the rotational symmetry of the original fleld 
theory will be restored on the lattice and the two notions of distance will only differ by 
an overall factor. 
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It is worthwhile pointing out that the simple choices (87) for the action 
and (88) for the measure lead to a partition function of the form 

7(n n 1 Na(T) Na-2{T) ^^^^ 

^{9d,9d~2) - 1^-^^ gd 9d-2 ' (y^) 

where = — log Kd and 5^-2 = logKd-2- Eq. (90) shows that the partition 
function is the generating function for the number of triangulations (of fixed 
topology) with given numbers Nd and Nd-2 of d- and (d — 2)-simplices. We 
thus reach the surprising conclusion that even in dimension d > 2, quantum 
gravity can be formulated as a (relatively simple) counting problem. 

4.4. INCLUSION OF MATTER FIELDS 

The discretization of matter fields coupled to dynamical triangulations is 
achieved by standard lattice field-theoretical methods. The simplest types 
of matter fields one can study are either scalar fields or (Potts) spin fields, 
carrying a discrete space-time label. One may also combine several fields 
of this type. The matter fields can be located at the vertices of the trian- 
gulation or at the centres of the d-simplices (that is, at the dual vertices). 
The interactions are typically of the form of nearest-neighbour interactions, 
where the "nearest neighbours" are the vertices that are one link length (or 
one dual link length) away from the original vertex. We expect these two 
formulations to become equivalent in the scaling limit. Some typical exam- 
ples of matter actions are 

^/ = f E(^-.-. -i)' (91) 

m 

where the ai are a set of Potts spins, or the Gaussian action 

S, = J2{<j,i-cj>jf (92) 

{ij} 

for a massless scalar field (p. Note that we did not need to include a coupling 
constant in front of the action (92). The massless scalar field can always 
be rescaled by a factor, which can then be absorbed by a redefinition of 
the coupling constants of the geometric sector. The coupling of Ising spins 
(Potts spins fjj with i = 1,2) and Gaussian fields to the 2d Lorentzian 
gravity model proceeds in a manner identical to (91) and (92). 

In higher dimensions we will also consider the coupling to gauge fields. 
As usual these are associated with the one-dimensional edges of the trian- 
gulation, and there are again two alternative formulations, depending on 
whether the links or the dual links are used. The gauge field action is more 
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complicated and we postpone its discussion to Section 6. The inclusion of 
fermionic degrees of freedom on a random manifold remains an open prob- 
lem, particularly in higher dimensions. It requires the definition of a spin 
connection on a simplicial manifold (c.f. [40]). The problem was solved in 
d = 2 [41] . In this case one can prove that a system of Wilson fermions on 
a triangulated manifold can be "bosonized" and represented as a system of 
Ising spins on the manifold [42, 41]. 

The discretized path-integral measure contains also a sum over matter 
fields. In the case of spin variables, the sum is simply taken over all pos- 
sible spin configurations. For continuous fields like the scalar fields above, 
one may consider non-trivial integration measures which introduce an addi- 
tional coupling to geometry. This possibility does not exist in field theories 
on fixed, regular lattices, where such a dependence is always trivial. It leads 
to some subtleties in the case of gauge fields, as we will discuss later. 

The simplest and most extensively studied example of a dynamically 
triangulated theory is that of Euclidean gravity on a two-sphere. The fun- 
damental building blocks in this case are equilateral triangles (2-simplices) . 
Triangles are glued together pairwise along edges (1-simplices), and each 
triangle has exactly three neighbours. The dual lattice is thus equivalent to 
a planar 0^-diagram. The curvature is localized at the vertices (0-simplices) 
which in general are shared by many triangles, each contributing 7r/3 to the 
total angle around the vertex. The regularity requirement introduced above 
prohibits configurations where a vertex is its own neighbour or where two 
vertices are connected by more than one link (in other words, closed loops 
of link length one or two are forbidden). 

The Lorentzian gravity model introduced in Section 2.1 may be viewed 
as a restricted version of the dynamically triangulated Euclidean model, 
since the triangulations contributing to the Lorentzian state sum constitute 
a subset of those appearing in the Euclidean system defined above. Recall 
also our construction of Euclidean from Lorentzian gravity in Section 3.1, 
by allowing for additional baby-universe branching. Again the set of all 
such geometries is a subset of all 2d simplicial manifolds, but both contin- 
uum theories coincide. What is at work in this latter case is "universality" , 
which ensures that the continuum limit is to a large extent independent of 
the short-distance details characterizing the class of triangulations we sum 
over. The universality properties of Euclidean 2d quantum gravity are well 
studied. For example, one may relax the manifold regularity condition to 
obtain a much larger class of simplicial complexes, whose continuum limit 
is still Euclidean quantum gravity. Only drastic modifications, like the sup- 
pression of baby universes, can move the system to a different universality 
class with a different critical behaviour and therefore a different continuum 
limit. Also around the fixed point leading to Lorentzian 2d gravity one 
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finds an independence of short-distance details. Universality with respect 
to a change of fundamental building blocks and the inclusion of higher cur- 
vature has been demonstrated in [28]. It is also encountered in a recently 
developped procedure for obtaining Lorentzian from Euclidean quantum 
gravity by removing baby universes [17]. There one ends up with a gener- 
alized class of triangulations (compared to the original Lorentzian model), 
but the continuum limit is still the same. 

5. Numerical setup 

5.1. MONTE CARLO METHOD AND ERGODIC MOVES 

Even in two dimensions, there are a number of issues of the matter-coupled 
theory that presently can only be addressed by numerical methods. Lorentz- 
ian gravity coupled to Ising spins has not yet been solved analytically. In 
matter-coupled Euclidean 2d quantum gravity, analytical considerations 
have not yet led to a determination of the fractal dimension of space-time 
(there are various suggestions leading to different answers), nor do we know 
what happens beyond the infamous c = 1 barrier, where analytical calcu- 
lations break down. In higher dimensions, we do not even have analytical 
solutions of the pure-gravity models. In all of these situations, numeri- 
cal simulations of the systems come in handy. They can answer specific 
questions and lead to unexpected results which in turn can inspire further 
analytical work. 

Numerical simulations of simplicial gravity have been the subject of a 
number of reviews (see [9, 43] for annual updates and [44] for more infor- 
mation on the computer codes used in the programs). Here we will only 
sketch the methods and use the simplest case of 2d Euclidean gravity as an 
illustration. Most of what we will have to say carries over to 2d Lorentzian 
gravity with only minor modifications. 

As explained in the previous section, the discretized theory can be de- 
scribed by the partition function 



where Seh is the discretized Einstein-Hilbert action (87), and the first 
sum is taken over all labelled triangulations of fixed spherical topology. For 
d = 2, the geometric part of the action simplifies and takes the form 



up to an additive constant proportional to the Euler number x, with fi = 
{k2 — kq)/2. In numerical simulations, it is simpler to use the labelled in- 
stead of the unlabelled triangulations. The partition function (93) is the 
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(93) 



Seh = IJ'N2, 



(94) 
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analogue of the grand canonical partition function in statistical mechanics. 
The cosmological constant Kd plays the role of a thermodynamic potential 
for the number of simplices. For general d, we may rewrite (93) as 

Z = J2e-'^-''^Z^„ (95) 

where ^jv^ can be interpreted as the partition function at fixed volume. 
A gravity-matter configuration C is uniquely specified by a geometry in 
the form of a labelled triangulation and by the values of all matter field 
variables. Each configuration enters the statistical sum with a probability 
proportional to e"'^^'-^), where S{C) = Seh + Smatter- As usual in statistical 
mechanics, physical information can be obtained by measuring the averages 
of various operators in this ensemble. The fact that each configuration has 
a real positive weight makes it possible to study the system (93) by Monte 
Carlo methods. The goal is to construct a numerical "generator" which 
produces configurations C with a probability -P(C) ~ e"'^^'-^). 

Except for very few cases, where a direct generation of all configurations 
is possible (e.g. for a conformal charge c = —2 in 2d), the standard way of 
obtaining such a generator is by means of a stochastic process (a Markov 
chain), which can be regarded as a random walk in configuration space. 
Since the random walk takes place on a computer, each step corresponds to 
the real time needed to perform such an operation. The stochastic process 
is defined by a function W{C — > C) > 0, giving the probability for a 
transition from a configuration C to C in one step. Assuming for simplicity 
that the configuration space is discrete, we have a normalization condition 

^ W{C C) = 1. (96) 
C" 

In most cases W{C ^ C) is chosen to vanish outside some "neighbour- 
hood" of C. Starting from an initial state pq{C), the system after n steps 
is characterized by a probability distribution pn{C), where 

Pn+iiC) = Y,PniC')WiC' ^ C). (97) 

C" 

The transition probability W must satisfy two basic requirements, 

1. ergodicity: any two configurations can be joined by a finite number of 
steps, and 

2. detailed balance: the condition 

P{C)W{C ^ C) = P{C')W{C' ^ C), (98) 
relating W{C C) to W{C' C). 
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It follows from (98) that W{C C) and W{C' C) are either both zero 
or both non-zero in which case they satisfy 

w{c C) _ pjC) 

W{C' C) ~ P{C) 

These two requirements guarantee that the stochastic process has a unique 
asymptotic probability distribution jJoo(C) ~ P{C), which is the only eigen- 
state of the transition matrix W with eigenvalue 1, 

J2poo{C')W{C' ^ C) = Poo(C). (100) 
c 

All other eigenvalues are strictly smaller than 1. This implies that - inde- 
pendent of the initial configuration - the system will reach the asymptotic 
distribution after infinitely many steps. Note that the asymptotic distribu- 
tion has the desired probability distribution. The rate at which the system 
approaches this limiting distribution depends on the other eigenvalues of 
W. The contributions from other eigenstates decay exponentially with the 
number of steps. The second-largest eigenvalue Ai provides us with an es- 
timate of the autocorrelation time T = — 1/logAi. When the number n of 
steps is 3> T, we can assume that the distribution is asymptotic. Typically 
the autocorrelation time behaves like T , where / counts the number of 
degrees of freedom of the system and (5 is a dynamical exponent, depending 
on the details of the algorithm. 

In a practical implementation the system starts from some configuration 
Cq. During the first step, it changes to Ci with probability W{Co — Ci), 
or remains at Co if the change is not accepted. After T steps it reaches 
a configuration C„ = C^^^ with a probability proportional to P{C^^^). This 
configuration is the starting point of a new process, during which another 
sufficiently large number of steps is performed and a new configuration C*^^) 
generated. Repeating this process we create a (finite) set of configurations 
{C^^\ . . . ,C^^^}, where depends on the computer time we spend on 
the project. The average of any operator in the ensemble of configurations 
defined by the probability distribution P{C) is approximated by an average 
over this finite sample of "typical" configurations. 

The requirements listed above by no means define the stochastic process 
uniquely. We are interested in efficient algorithms which minimize 8, and 
can produce a large number of independent configurations in the shortest 
possible time. There is no simple way to guess at the outset whether an 
algorithm is efficient or not. Each problem must be treated individually 
and the autocorrelation time measured numerically. There are some gen- 
eral guidelines which one usually follows when creating a new algorithm. 
As discussed above, the efficiency depends on the choice of the transition 
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probabilities W(C ^ C). We would like the algorithm to have a high 
"mobility", that is, a high probability that a configuration of the system 
will change at each step. This means that the set {C} of configurations 
which can be reached from a given C should be limited. If it is too large, 
each transition probability will be small, implying that the configuration 
will most likely not change. On the other hand, the set {C} must be large 
enough to ensure ergodicity. 

The detailed-balance condition (99) implies that if both transitions are 
to be reasonably probable, we must choose the set {C} such that its el- 
ements have similar probabilities to that of C (equivalently, have a small 
action difference \S{C) — S{C')\). On a fixed lattice, small differences in the 
action are usually realized by considering at each step only local changes in 
the field variables, for instance, changing only one variable, while keeping 
all the others fixed. When a more complicated change is attempted, the 
action difference is in general large and proportional to the volume of the 
system. Local changes are not very efficient when the typical fluctuations 
are long-ranged, as happens close to a continuous phase transition. Creat- 
ing a Monte Carlo algorithm which is at the same time ergodic and has a 
reasonably small autocorrelation time even in the critical region is an art. 

The first step in constructing an algorithm for simplicial gravity is to de- 
fine a method of coding the configurations. Prom a numerical point of view 
it is natural to work with labelled rather than unlabelled triangulations, 
because otherwise it is almost impossible to keep track of the (dynamical) 
connectivity. To code the geometric structure of such a configuration it is 
in principle sufficient to have a list of the vertex labels of all simplices (of 
all dimensions) of the triangulation. Two vertices are neighbours if they 
belong to the same simplex. From this list we can reconstruct the complete 
geometry of the system. The fact that we have a list of simplices means in 
practice that we have also assigned labels to the simplices. Two simplices 
are neighbours if they share a {d — l)-simplex or, equivalently, d vertices. In 
addition, we can make a list of all subsimplices and count their order. The 
process of reconstruction may still be complicated and it is often useful to 
keep at each step even more information, in the extreme case the lists of 
all subsimplices. The more information we keep, the easier it is to recon- 
struct local properties of the geometry. However, it also means that during 
geometry updates more data will have to be changed. 

The next step consists in constructing the stochastic process described 
above. In our case this amounts to defining a set of "moves" , which connect 
configurations with different geometry and matter content. Usually these 
two changes are performed separately: changes of the matter fields are per- 
formed for a fixed geometric configuration, by techniques very similar to 
those used for fixed lattices. For spin systems there exist highly effective 
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cluster algorithms [45] with very short autocorrelation times. Cluster algo- 
rithms are special in that they permit global changes of the spin configura- 
tions. For Gaussian fields, no such algorithms exist and the field variables 
are usually updated sequentially. (There are some methods to shorten their 
autocorrelation time.) Generally speaking, the autocorrelation times for 
Gaussian fields are much longer than those for spin systems. It is therefore 
necessary to repeat the updating of all fields on the manifold several times 
in order to obtain independent field configurations. 

In addition, we must define a set of geometric moves which generate 
local changes of the geometry, while preserving its underlying topology. 
Such moves are necessary even in the absence of matter. "Local changes" 
means that only a small region of the manifold is affected at each step. 
Again, it is important that the set of moves be ergodic in the space of 
all possible triangulations. For simplicial manifolds, one possible choice is 
given by the so-called Alexander moves [46] . In practical applications most 
algorithms for a d-dimensional geometry are based on a finite subset of 
moves containing d + 1 operations. We will describe this set of operations 
in the case of 2d and postpone its generalization to higher d to Section 6. 

The first operation is called a flip and involves two neighbouring tri- 
angles. The triangles are denoted by their vertex labels, {123} and {124}, 
and share the link {12}. It is important that the four labels are all distinct 
(excluding {3} = {4}, say) and that the vertices {3} and {4} are not con- 
nected by a link. The flip consists in replacing this configuration by the 
two triangles {341} and {342}. In other words, the link {12} is "fiipped" 
to the link {34}. The restrictions imposed above guarantee that this move 
does not produce a pathological triangulation. The inverse of the flip oper- 
ation is again a flip. Note that the flip move takes place entirely inside the 
link loop {1324}, whose boundary it leaves unchanged. Since two simplices 
are replaced by two others, the fiip is sometimes called a (2, 2)-move. A 
fiip can be performed almost everywhere in the manifold, giving rise to a 
different labelled manifold whose connectivity is changed locally. If matter 
fields are present at the vertices, one usually assumes that their values are 
unchanged. The transition weights W{C C) and W{C' C) for a flip 
and its inverse can easily be computed from the detailed-balance condition. 
We will not be more specific here, since this depends on the details of the 
geometric coding (in particular, on how links are selected). A flip move 
leaves the numbers of triangles and vertices unchanged. 

The second move adds a new vertex {4} inside a triangle {123}. This 
move replaces the old triangle by the three triangles {124}, {234} and 
{314}, and leaves the manifold outside the "boundary" of {123} unaffected. 
It is also known as a (1, 3)-move. The new vertex is of order three, since 
it is shared by three triangles. The third move is the inverse of the (1,3)- 
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Figure 4- A sot of three moves which is ergodic in the class of two-dimensional trian- 
gulations of fixed topology. The first diagram shows the (3, l)-move and its inverse. The 
second diagram shows the (2, 2)- or flip move. By itself, this move is ergodic in the class 
of triangulations of fixed volume N2 and topology. 

move. It is a (3, l)-move where a vertex of order three is removed, and 
the three triangles which share it are replaced by a single triangle. The 
(1, 3)-move can be performed on each triangle of a manifold, but its inverse 
needs a vertex of coordination number three. It is obvious that in both 
cases the configurations before and after the move are two different labelled 
triangulations. If matter fields are present, a new field must be created at 
the new vertex generated during a (l,3)-move. This must be taken into 
account when the transition weights are calculated. The three types of 
moves are depicted in Fig. 4. 

The three operations just described leave the manifold topology un- 
changed. They form an ergodic set, which means that any two simplicial 
manifolds of the same topology can be related by a finite sequence of moves. 
In 2d one can construct such a sequence explicitly. The set of ergodic moves 
is not unique, and alternative sets of local moves have been used in appli- 
cations. For example, a point-splitting algorithm is described in [8]. The 
point-splitting move and its inverse (illustrated in Fig. 5) change the volume 
N2 by ±2. The (1, 3)- and (3, l)-moves are special cases and the (2, 2)-move 
can be realized as a sequence of two point-splitting moves. 

Our discussion so far suggests that we may set up a numerical simulation 
by generating configurations according to the probability distribution (93) , 
and use the resulting sample of configurations to measure the quantities of 
interest. However, this is not really feasible, since the system described by 
(93) is an open system in the sense that arbitrarily large configurations may 
be produced by using one of the sets of geodesic moves described above. 
In practice, we must limit this size because of the obvious memory restric- 
tions of a computer. A simple solution is to generate a set of configurations 
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Figure 5. The point-splitting moves constitute an alternative set of ergodic moves. 




Figure 6. The move used in the Monte Carlo updating of a 2d Lorentzian geometry. 

of fixed volume N2, and repeat the experiment for various values of N2. 
This is what one typically does in numerical simulations of field theory. For 
the case of Euclidean 2d triangulations we are particularly lucky, since the 
(2,2)-fiip move is already by itself ergodic in the space of all triangulations 
of fixed volume N2, considerably simplifying the computer simulations. Un- 
fortunately there is no similar result in higher dimensions; we will describe 
later the method used in this case. The flip move cannot be used in numer- 
ical simulations of the Lorentzian model described in Section 2.1, since it is 
not compatible with the causal structure. One uses instead a version of the 
point-splitting move (Fig. 6) to update the geometry (see [12] for details). 

We conclude this section by describing another type of move which 
is associated with a large change in matter and geometry (while keeping 
AS* small). It is motivated by the baby-universe structure of Euclidean 
simplicial gravity. Let us consider a two-dimensional triangulated manifold 
with the topology of a disc. It is characterized by its volume n and boundary 
length I. For a random geometry these two numbers are not related since 
for any finite I, n can become arbitrarily large. If we have two such discs, of 
volume n and N — n, but with identical boundary lengths, we can identify 
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their boundaries to create a closed simplicial complex of spherical topology 
and volume N2 = N. The smaller component is the baby universe and the 
larger one the parent universe. A typical spherical manifold observed in 
numerical simulations will contain many structures of this kind, of various 
sizes and lengths of boundaries, called "necks" . The shortest possible neck 
has length three, and the associated baby universe will be called a minimal 
baby universe (minbu). We may also have baby universes characterized by 
longer (but finite and small) neck sizes. 

The existence of baby universes opens a completely new possibility for 
updating the matter sector of the theory. Note that for a finite boundary 
length I, the interaction between the matter in the baby universe and the 
parent universe will be of order I. Take the example of Ising spins: even 
if we flip all spins inside the baby universe, the change in the action will 
only depend on the interactions across the boundary, and not the size n of 
the baby universe. This is a completely different situation than on a regu- 
lar lattice where large domains always have large boundaries. We may use 
this observation to define a matter update which induces a large change in 
magnetization and has a large acceptance. The efficiency of such an algo- 
rithm will depend on the typical baby universe size. One possible algorithm 
consists in searching the triangulation for a minimal neck and flip all spins 
inside the associated minbu with probability W. The same technique can 
be applied to massless scalar fields. Here we can perform two operations 
on the fields in the minbu: change all of their signs and/or add the same 
constant to all of them. 

The existence of minbus can be exploited also in the geometric sector. 
For example, we may use the following three-step algorithm. 

• Locate a minimal neck in the triangulation, and cut open the original 
spherical manifold along the neck. Both of the resulting discs are of 
the form of a triangulated sphere with a triangle removed. 

• Close off both holes by a single triangle, to obtain two spherical trian- 
gulations of volumes n + 1 and N — n + 1 (if the original triangulation 
had volume N). 

• Remove an arbitrary triangle from each of the two manifolds, and glue 
the resulting discs along their triangular boundaries. 

We must compute the correct probability factors for such a move and 
eventually also include changes in the matter fields. Based on this idea, 
an extremely efficient "minbu surgery" algorithm has been constructed, 
shortening the autocorrelation times by three orders of magnitude [16]. The 
move does not even require a lot of geometric updating if it is accepted, 
since the connectivities of the manifold are changed only locally. The minbu 
surgery usually supplements the local moves described earlier, although in 
2d one may construct an algorithm that is based exclusively on the cutting 
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and pasting of baby universes, if one permits also longer neck lengths. (A 
flip move may be viewed as a particular realization of this, if instead of 
minbus we consider baby universes with boundaries of length I = 4.) 

5.2. OBSERVABLES IN 2D EUCLIDEAN GRAVITY 

The Monte Carlo method described above can be used to generate a set of 
uncorrected configurations of volume N2 which are distributed according to 
their Boltzmann weights . In order to understand the finite-size effects, 
one needs to perform the simulations at different values of N2- The sizes 
that have been used in numerical simulations to date range from 1000 to 
128000 simplices. 

In 2d it speeds up the calculations to extend the class of allowed trian- 
gulations to include also loops of link length one and two (corresponding to 
self-energy and tadpole subdiagrams on the dual lattice). Let us denote this 
extended class by T and the corresponding class of labelled triangulations 
by T. Although the complexes constructed in this way are no longer sim- 
plicial manifolds, we can still keep the notions of global topology, of local 
neighbourhoods and of a geodesic distance. Likewise, the relations between 
the numbers Nj. remain unchanged. In some cases (including pure 2d grav- 
ity), models based on this set of geometries can be compared directly with 
the analytic solutions of corresponding matrix models, where the exclusion 
of tadpole and self-energy subdiagrams corresponds merely to a finite renor- 
malization of the bare coupling constants. On T, the numerical simulations 
become even simpler. For fixed volume, the flip move is still ergodic and 
also the minbu surgery moves can be generalized. The advantage of using 
this class of triangulations is a reduction of the flnite-size effects, since it 
turns out that the local restrictions on the connectivity do not affect the 
scaling properties of the system. 

Our next step will be to describe the measurement of suitable "observ- 
ables" on the ensemble of configurations generated by the Monte Carlo 
algorithm. The observables most easily obtained are the critical exponents 
related to the geometry or the matter fields. In two dimensions, one can 
sometimes obtain such observables analytically, and use them to test the va- 
lidity of the numerical results. Recall that in 2d we start from the partition 
function 



where Zy is the partition function at fixed volume. If the central charge of 
matter is c < 1, it can be shown analytically that Zy behaves like 



f <t>i V 



(101) 



Zv ~ e^°^y^=*'-3(l + 0(1/F)) 



(102) 
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for large V and spherical topology. The subleading power contains the crit- 
ical exponent ^ystr, which has a known dependence on c (see formula (125) 
below). The pure-gravity case corresponds to 7^^^ = ~l/2- 

There are other statistical systems whose partition function behaves 
like (102), most notably, various realizations of branched polymers [20, 47]. 
In these models, jstr is positive and < 1/2, and V counts the number of 
vertices. 

Measuring the distribution of minbu sizes for a triangulation of fixed 
volume N2 = y is an efficient device for determining the critical exponent 
jstr [48, 49]. The average number {b{n))y of minbus with volume ^ n < 
V/2 is given by 

{b{n))y — , (103) 

since a minbu of size n can be regarded as a spherical triangulation of 
volume n + 1 with one marked triangle. (Since both n and V — n are 
assumed large we can neglect small corrections to the volume.) Using (102), 
we obtain 

{h{n))y ~ V {n{V - n))^^*^-^ . (104) 

Measuring ^str gives us information about the fractal structure of the 
theory and provides a simple test of algorithms in cases where its value is 
known. In all cases with c < 1, one finds an excellent numerical agreement 
with the predicted values. Prom relation (104) wc sec that for —1 < ^str < 0, 
the average minbu size {n.)y remains finite, but that {n^)y ~ +1. This 
means that with growing ^gtr we will observe increasingly large fiuctuations 
in minbu size. If jstr > the average minbu size behaves like {n)y ~ V^"*'^ . 

Another important observable that can be studied by numerical simu- 
lations is the volume-volume correlator, which is a particular example of a 
geometric two-point function. In the continuum theory, it can be defined 
as 

Ga{R) = J V[g^,]V<j> e-^[^''^l j I d^de' V5(0 S{dg{^, O - R), 

(105) 

where dg{^,,^') denotes the geodesic distance between two points ^ and 
calculated with respect to the metric 5^1/ . Ga{R) is the partition function 
for universes with two marked points separated by a geodesic distance R, 
which we already encountered in Section 3.2. In a discretized theory we can 
construct the analogous quantity 

G^r) = E E e-'^^^-^— (^^) E SD(i,j),r, (106) 
r <t>i hi 

where D{i,j) is now one of the integer-valued geodesic distances introduced 
at the end of Section 4.2. The two possibilities will in general differ by a 
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finite scaling factor. For large r, one can show that Gn{r) falls off exponen- 
tially like 

G^(r) ~ e-"*('')^ r » l/m(/i). (107) 
For small r we expect 

G^(r) ~ (108) 

where ?7 is the anomalous dimension of the two-point function. For fi ^ /iq 
we make the ansatz m(/x) ~ (/i — fxo)^^'^'^ where dh is another critical expo- 
nent of the theory. Prom the functional form of the partition function G^{r), 
one can derive the estimates {V) ^ ~ — /xq) and (r)^ ~ — fio)^/^^. 
From these one obtains a relation between the average linear extension and 
the average volume of the configuration, namely, 

(V)^ - {rf; . (109) 

The exponent dh is called the cosmological HausdorfF dimension. It is a 

large-scale property of the average "quantum geometry" of the ensemble, 
and therefore need not coincide with the dimension of the individual trian- 
gulations. 

The function Gn{r) has been calculated analytically for pure gravity, 
as explained in Section 3.2, leading to d/j = 4. As we have already pointed 
out, it is convenient in numerical simulations to work with triangulations 
of fixed total volume. The two-point functions at fixed and unrestricted 
volumes are related by a discrete Laplace transform. From definition (106) 
we have 

G'M(0 = Ee"^''G'y(r), (110) 
V 

where Gv{r) is the partition function at fixed volume V. For a particular 
configuration of the system, Gy{r) is measured as follows. Start with a 
(dual) vertex i and find all (dual) vertices at (dual) distance 1 from it. By 
iterating this process, the triangulation is decomposed into shells character- 
ized by their distance r from the initial point i. Note that the shells will in 
general be disconnected. We can measure the total length of each shell, that 
is, the number of vertices in the shell.^ This construction is repeated for all 
starting points and all configurations in the sample to obtain an averaged 
distribution. The resulting quantity (n(r))y is - up to a normalization - 
a numerical estimate of Gv{r). It is convenient to use the normalization 
condition 

J2Gv{r) = V, (111) 

r 

®Note that one could also measure the number of disconnected parts and/or the higher 
moments of their length distribution. 



QUANTUM GRAVITY 



41 



which leads to an interpretation of Gy(r) as the average volume of a shell 
with radius r for a triangulation of volume V. For large V one expects that 

Gv{r) ~ V^-^/'^'^F{x), (112) 

where x = r jV'^l'^^ is the scaling variable. In terms of x, we have a normal- 
ization condition 

\xF{x) = \. (113) 



Prom an inverse Laplace transform of relation (107), we can deduce that 
for large x, F(x) must behave like 



logF(x) ~ x^h^ x»l. (114) 

On the other hand, wc expect for small x that 

F{x) - (115) 

where du is another Hausdorff dimension characterizing the quantum ge- 
ometry at short distances, and is related to the anomalous dimension ry. 
Measurements of Gy(r) have been performed for systems with a variety of 
matter types. A simple way to determine the cosmological Hausdorff dimen- 
sion dh is to plot Gy(r) against r. Even if the function F{x) is not known 
explicitly, one can obtain an estimate of dh by comparing the maxima of 
these curves for different volumes V . One then plots the rescaled quantities 
GY(r')lV^~^l'^^ as a function of a;, where according to (112) the curves for 
different volumes V should fall on top of each other. This procedure can 
be viewed as a finite-size scaling method, of the type used in the study of 
critical phenomena in statistical physics. In the context of quantum gravity, 
it was introduced in refs. [50, 51] for d = 2, and in rcfs. [52, 53] for d = A. 

For pure two-dimensional gravity one can compute F{x) exactly, and 
finds dff = dh = 4. Comparing this to the numerical analysis, it was discov- 
ered that a scaling form with the correct value of d^ can only be obtained 
by using a shifted scaling variable defined by 

r + 5 

= (116) 

where the constant 6 is determined numerically. In the case of pure gravity 
the shift 5 can be calculated analytically, but for more complicated systems 
this is usually not the case. The agreement between theory and numerical 
simulations is impressive. In Fig. 7 we show the scaling of Gy {r) for various 
volumes V and the theoretical curve (which can hardly be seen since it 
coincides perfectly with the numerical data). 
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Figure 7. The distributions n{x) = Gv{x)/V^ ^^'^'^ in 2d Euclidean gravity, for volumes 
V = 1000, 2000, 4000, 8000, 16000 and 32000, together with the theoretical distribution 
F{x). 

Such a shift is also necessary in more compHcated systems. It was soon 
realized that for the case of pure gravity this shift is a leading-order finite- 
size correction to the scaling form (112), which was obtained for V oo. 
This improved scaling could also be used to check the small-x behaviour of 
F(x) with very good accuracy. 

Another important test of our numerical methods is to check whether 
the scaling limit of the two-point function is indeed independent of the de- 
tailed definition of the distance function. Numerically it is easy to measure 
the two-point function Gv{r), where r is now the link distance on the tri- 
angulation. In this case we have no exact analytical prediction, but we may 
again fit the numerical results (which look very similar to those for Gy(r)) 
to the scaling relation (112). After a trivial rescaling of the distance, the 
matching with the theoretical curve F{x) is even better than previously, 
but again one must include a shift 5. Numerically, the shift is much smaller 
than before. This was to be expected, since a shift in x compensates for 
small-distance artifacts, which are larger on a lattice with smaller vertex 
valence. The aspects mentioned above are illustrated in Fig. 8 by data 
obtained in the study of a c = —2 conformal field theory coupled to 2d 
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Figure 8. The two upper figures depiet the measured distributions Gv ij') of 2d Euclidean 
gravity coupled to c = — 2 conforrrial matter, defined in terms of the dual link distance 
and rescaled according to eq. (112) to get F{x), with a shift 5 = 4.5 (left) and without a 
shift (right). In the two lower figures, the link distance has been used, again with a shift 
5 = 0.5 (left) and without (right). The best fit for the Hausdorff dimension dn extracted 
from the data is 3.58 ± 0.01, in agreement with theoretical predictions [55]. 



Euclidean gravity [54]. 

Analogous rcstilts hold for branchcd-polymer systems, whose two-point 
functions can be calculated exactly. The distance r is in this case given 
by the (unique) number of links separating each pair of points. For fixed 
volume V — ^ oo, we find a universal scaling of the form (112), but now with 
an (intrinsic) Hausdorff dimension = 2, which also coincides with the 
value for ■ Again a shift in x is needed as a leading finite-size correction 
to obtain optimal agreement with the theoretical results. 

To summarize, measuring Gy(r) or (5y(r) provides an efficient way for 
determining the Hausdorff dimension dh, also in cases where the two-point 
functions are not known exactly. To take care of short-distance effects, it is 
important to work with an improved, shifted scaling variable xs. In practice 
one measures Gy(r) for a number of volumes and tries to find a best fit 
to the universal scaling relation (112) by fitting both 5 and dh- 

Another geometric observable characterizing a compact metric manifold 
is its spectral dimension, which is related to the diffusion equation and the 
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spectrum of the Laplace operator. If a diffusion process is started with a 
completely localized initial condition, one can measure the return probabil- 
ity P{T) of a fictitious test particle after "time" T. On a fixed continuum 
geometry, one obtains after averaging over the initial point that 

P(T)-^J:^'^T\ (117) 

k 

for small T, where the coefficients can be expressed in terms of local 
geometric invariants. We can measure dg also in the discretized theory, but 
some care must be exercised when comparing with a continuum formula like 
(117), since it is exactly the short-time limit T — that is ill-defined in the 
discrete case. This can already be seen in the simple case of diffusion on a 
discrete one-dimcnsional line, with a discrete Laplacian and a continuous 
time. The diffusion equation 

#i (pi+i + (f)i-i - 2(f)i 

df 2 ' ^^^^^ 

with the initial condition ^i(O) = 6io can be solved exactly, yielding 

P{T) = e-^Io(r) - ^(1 + 0{l/T)), (119) 

where Io(T) is the Bcsscl function. One therefore rederives the correct spec- 
tral dimension (dg = 1) only for large T. The short-time behaviour is com- 
pletely different and dominated by discretization artifacts. For this simple 
example, one can check numerically the effect of the finite system size, by 
changing the discrete line to a closed circle with V points. One observes 
three distinct regions: 

• small time T, which is similar to the example above, and dominated 
by discretization effects; 

• intermediate T, where one obtains the correct value of dg; and 

• very large T, where the system approaches a stationary translation- 
invariant state. 

A similar structure is also expected in less trivial cases, such as simplicial 
quantum gravity in two dimensions and higher. This behaviour is illustrated 
in Fig. 9, where we show the results of measuring the spectral dimension for 
2d Euclidean quantum gravity coupled to matter fields of various central 
charges c. 

The measurement of the spectral dimension in 2d Euclidean quantum 
gravity is a nice illustration of the fruitful interaction between numerical 
"experiment" and theory. Since the quantum geometry of 2d Euclidean 
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Figure 9. The spectral dimension ds ~ — 2cilog P(r)/dlog T versus T for c = (top 
curve), c = 1/2, c = 1, c = 3 and c=5 (bottom curve) theories coupled to 2d Euclidean 
quantum gravity. The system size is ^=16000 triangles. 

gravity has a fractal HausdorfF dimension = 4, one eould also have ex- 
pected an anomalous spectral dimension dg. However, numerical simula- 
tions consistently found that = 2, no matter which field theory was 
coupled to gravity, as long as the central charge c remained < 1. This in- 
spired an analytic proof of this value, using continuum conformal Liouville 
theory and an assumption about finite-size scaling [56, 57]. We should point 
out that ds = 2 does not imply that one cannot see the fractal structure 
corresponding to d/j = 4 in the diffusion process. It manifests itself in the 
relation 

{R)t ~ T^/'^'' (120) 

between the "time" T and the average geodesic distance R diffused in that 
time. We refer again to [56, 57] for details. 

For c > 1, the numerical estimate was close to dg = 4/3. Eventually it 
was proven that generic branched polymers have a spectral dimension ds = 
4/3 [58], thus reconciling theory and numerical results. At the same time, 
it provided independent evidence that for c > 1 the space-time degenerates 
into branched polymers. 

We have so far discussed only purely geometric observables. For theories 
coupled to matter, there will be additional critical exponents characteriz- 
ing the behaviour of the system close to the critical matter coupling f5(.- 
Near this transition point, we expect to find a divergent correlation length 
^ \I3 — Pc\~'', for some positive u. Since we always study systems of 
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finite size, we will never observe a genuine phase transition, but merely a 
pseudo-critical point /?* (V) , where the linear extension L of the system be- 
haves like L ~ For a finite system L*^ = V. A system with a dynamical 
geometry may possess two scales: one associated with the cosmological con- 
stant when /i ^ /io and another one related to the matter phase transition. 
Using standard finite-size scaling arguments, we expect that the measured 
values of the specific heat c^, the magnetization per unit volume m, and 
the magnetic susceptibility x at the pseudo-critical point /?* behave like 

mm ~ 

Xm ~ (121) 

where a, (3^ and satisfy the scaling relations 

a + 2l3m + lm = 2, 2pm + 7m = T^d. (122) 

For magnetic systems corresponding to < c < 1, one knows from analyti- 
cal calculations the values of the critical exponents as well as the values of 
the critical points fic- 

A relation between the scalings in the geometric and matter sectors can 
be obtained by generalizing the two-point functions G^{r) and Gy(r), to 
include also correlations between the matter fields at geodesic distance r. 
One example of a quantity of this kind is 

which may be expressed as the discrete Laplace transform of a finite- volume 
correlator, 

Gl^'^^r) = Y^e-^'^Glf'^Hr). (124) 

V 

In this expression, 0j and ipj are local functions of the matter or geometric 
variables at the points i and j. A typical example is the spin-spin correla- 
tion function of the Ising system. At the combined critical point of the Ising 
model and the cosmological constant, (123) yields a reparameterization- 
invariant definition of a matter correlator in quantum gravity. More pre- 
cisely, the correlation function is that of a c= 1/2 conformal field theory 
coupled to 2d Euclidean quantum gravity. Note also that the previously 
defined geometric two-point function (106) is a special case of (123), with 
the unit operator replacing the spin operators. 

This construction has provided us with a new understanding of the 
KPZ-exponents of conformal field theories coupled to gravity. These ex- 
ponents may be viewed as the "dressed" scaling exponents of the primary 
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Figure 10. The normalized spin-spin correlation function for the Ising model coupled 
to 2d Euclidean gravity, plotted as a function of the dimensionless length x = r/N^'^'^ , 
where r is the geodesic length and Nt the space-time volume. The number of triangles Nt 
ranges from 2000 to 256000 and the Hausdorff dimension is dh = 4.0. The short-distance 
behaviour is in accordance with KPZ scaling. 



fields in the conformal theory. However, before it was reahzed that they 
should be analyzed in terms of the definition (123), it was not clear how 
to generalize the usual fiat-space correlators because of the requirement 
of reparamctrization-invariance (see [59, 60] for further discussion). Fig. 
10 shows the result of a numerical simulation of the (suitably normalized) 
spin-spin correlator. 

We should point out a subtlety in the definition of the correlator (123). 
If the expectation values of the fields entering are not zero, one might 
be interested in defining the corresponding connected correlator. However, 
there is no straightforward way of defining such an object since (123), in 
addition to the usual average over field configurations, also includes an 
average over geometry. A discussion of possible definitions of the connected 
part of the two-point correlators and their scaling properties can be found 
in [61]. 

All of the observables introduced above can also be used in numerical 
simulations of higher-dimensional gravity-matter systems, and - apart from 
the ones involving baby-universe counting - also for Lorentzian quantum 
gravity. 
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5.3. COMMENTS ON THE 2D RESULTS 

Our obvious starting point in 2d Euclidean quantum gravity was the mea- 
surement of ^str , in order to compare it with the theoretical prediction for 
c < 1, namely, 

c-1- V(c-l)(c-25) 

Istr = J2 • ^'-^^^ 

For c > 1, this quantity becomes imaginary and its interpretation in terms 
of Liouville theory breaks down. For c = 1, logarithmic corrections appear. 
For a number of special values of c, we have explicit discrete models, which 
at their critical points represent conformal field theories with charge c. They 
are pure 2d gravity (c = 0, ^str = the Ising model (c = 1/2, jstr = 

— 1/3), the 3-state Potts model (c = 4/5, jstr = —1/5), the 4-state Potts 
model (c = 1, jstr = 0), cind a single massless scalar field (c = 1, jstr = 0)- 
In all of these cases, formula (125) has been verified with great accuracy. A 
coupling of several matter fields corresponds to adding their central charges. 

Away from their critical points /3c, the finite spin systems approach a 
pure-gravity behaviour, with jstr —1/2. This observation can be turned 
into a method for locating /3c, by monitoring the change in jstr ^ the 
couphng /3 is varied (Fig. 11). 

For all the spin systems, measurements of the magnetic critical ex- 
ponents defined in (121) arc in perfect agreement with the theoretical 
predictions^. We conclude that for the matter models with < c < 1 
which have been used in numerical simulations, the scaling limit of the dis- 
cretized theory corresponds to a unitary matter model coupled to Liouville 
gravity. 

A somewhat surprising result in this range of c is the numerical value 
of dh obtained from the scaling of Gv{r). For < c < 1, there are two 
theoretical predictions for dh, namely 



d« = 2^^+^, (126) 

derived by studying diffusion in a fluctuating geometry, using Liouville 
theory [55], and 

^ (x/r^+x/25^)' ^ ' 

'^This also includes the case c = 1, which can be realized either as two Ising spins, as 
a g = 4 Potts model, or as a single Gaussian field coupled to gravity. However, there are 
logarithmic corrections to formula (102) which have to be included in the fits in order to 
obtain ^str = from the data. 
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Figure 11. Measurements of •ystr as a function of (3, for the Ising model (white), and 
the 3-state Potts model (black), at system sizes N = 1000 (circles) and 2000 (squares). 

derived from matrix model considerations [62]. The formulas agree for c=0 
(where = 4) , but differ elsewhere. In both derivations there are plausi- 
ble, but unjustified, assumptions. Neither of them could be confirmed by 
numerical simulations®. Instead, a numerical estimate d/^ ~ 4 was obtained 
for all cases with < c < 1 [51, 59] together with dn — dh [60, 59] from 
the small-x behaviour of F{x). The numerical errors seem sufficiently small 
to support the conjecture that dh = ^ for unitary matter theories coupled 
to 2d gravity. The spin-spin correlation functions measured for the Ising 
system (c = 1/2) at /? = /?c seem to indicate that the matter degrees of 
freedom do not introduce a new scale and that also for these correlations, 
there is a universal scaling corresponding to d^ ~ 4. 

Our numerical information for the case c < is only partial. A special 
case is c = —2, where the simulations can be based on a direct generation 
of diagrams, instead of Monte Carlo methods. As already mentioned in 
connection with Fig. 8 this led to dh = 3.58 it 0.02, which is very close to 

*It is now understood that formula (127) measures the fractal dimension of (general- 
ized) spin clusters of the matter fields, rather than the fractal dimension of the underlying 
geometry [63]. 
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the value of dh = 3.561 . . . predicted by (126). Less conclusive numerical 
data for c = — 5 also show agreement with (126). Formula (126) also predicts 
that dfi = 2 as c ^ — oo, which is the expected flat-space behaviour. It thus 
seems that eq. (126) could be correct for c < 0. It is still not understood 
why the formula fails for < c < 1. 

For c > 1, the minbu size diverges as F — oo. Systems with c > 1 can 
easily be studied numerically by coupling several copies of spins or scalar 
fields. In all cases studied, ^str was found to be positive and to approach 
+1/2 for c ~ 4 and bigger, suggestive of a branched-polymer system. This 
interpretation would also agree with the scaling analysis of G'y(r) which 
is consistent with dh ^ 2, and the measurement of dg, consistent with 

— >■ 4/3. It is not clear to what extent the observed behaviour is universal, 
since the numerical estimates depend on the details of the regularization. 
For example, the approach to the branched-polymer phase seems slower 
when regular triangulations T rather than generalized configurations in T 
are used. This could be due to finite-size effects, or alternatively to the 
presence of other universality classes of branched-polymer systems. 

In the 2d Euclidean case, we have learned from both theoretical and 
numerical investigations that there is only a small "window" of physically 
sensible theories, by which we mean unitary matter models coupled to grav- 
ity. If wc move oTitsidc this range by adding too many matter fields, the 
theory breaks down. This indicates that matter and geometry are interact- 
ing strongly. 

Baby universes play a crucial role in understanding the nature of the 
spin-gravity interaction, highlighting at the same time the difference be- 
tween the Euclidean and Lorentzian gravity models. Again the Ising model 
serves as an ideal illustration. The Ising ground state at zero temperature, 
for both fixed and fluctuating geometries, is the state where all spins are 
aligned. Since the energy of a given spin configuration is proportional to 
the length of the boundary separating spin-up and spin-down regions, the 
dominant spin configurations at low temperature (equivalently, at large spin 
coupling) are those with minimal spin boundary lengths. 

The lowest spin excitations contributing to the free energy density come 
from spin clusters with boundaries of minimal length. Unlike on flat, regular 
lattices, a short spin boundary does not imply that there are few spins 
inside. On the contrary, the typical situation in Euclidean quantum gravity 
has an entire baby universe of up-spins, say, on one side of the minimal 
boundary, and the parent universe with opposite spin orientation on the 
other, as illustrated in Fig. 12. Since the baby universe can have any size, 
there is no restriction on the size of spin clusters, even at low temperature. 
The introduction of a fluctuating geometry with baby universes thus has a 
strong effect on the matter behaviour. In turn, since short spin boundaries 
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Figure 12. Two spin clusters separated by only two links which form a "bottleneck" on 
the surface. Two of the four triangles containing one of the two spin-boundary links are 
also shown. 




Figure 13. A surface with a minimal boundary and its recursive decomposition into 
baby universes associated with spin clusters. 

are energetically preferred, the matter has a tendency to "squeeze off" the 
underlying geometry, generating even more baby universes. 

When more than two Ising models are coupled to 2d Euclidean quantum 
gravity, the matter-geometry interaction becomes so strong that the geome- 
tries degenerate into so-called branched polymers, which can be viewed as 
infinitely branched trees of baby universes of cut-off size (the lattice spac- 
ing a). This provides us with an explicit picture of the c = 1 barrier of 
2d Euclidean quantum gravity. The proliferation of baby universes and its 
relation to spin clusters arc illustrated in Fig. 13. 

Since baby universes are absent from Lorentzian gravity, its interaction 
with matter is much weaker. The Ising model on Lorentzian geometries has 
not yet been solved analytically, but its critical exponents have been deter- 
mined both by a high-temperature expansion and Monte Carlo techniques 
[12]. The Hausdorff dimension of space-time is still dh = 2, and the critical 
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Figure 14- Two typical space-time configurations, resulting fi-om the coupling of one 
Ising model (left) and eight Ising models (right) to Lorentzian gravity. 

Ising exponents retain the Onsager values found on fixed, regular lattices, in 
spite of large fluctuations of the geometry. It has also been shown [28] that 
a particular dimer model^ coupled to Lorentzian gravity does not change 
the fractal dimension of the geometry. It is thus tempting to conjecture 
that the (fractal) dimension of space-time remains equal to 2 as long as the 
central change of the matter fields is less than or equal to 1. 

It is an interesting question whether there is an analogue of the c = 1 bar- 
rier when sufficiently many matter degrees of freedom (with a sufficiently 
large central charge c) are added. This is indeed what seems to happen. A 
phase transition in the geometry has been observed when coupling 8 Ising 
spin copies (corresponding to c=4) to 2d Lorentzian quantum gravity [13]. 
This effect is illustrated in Fig. 14, where we show typical configurations of 
the fluctuating Lorentzian geometry in the case of coupling to a single Ising 
model and to 8 Ising models. The one-Ising triangulation is qualitatively 
very similar to the fluctuating geometry without any Ising spins, but in 
the case of 8 Ising spins, the effect of the matter is very pronounced. Since 
the creation of baby-universe branches is by definition forbidden, the only 

^It is known that the critical behaviour of dimer models is associated with a c = 
—2 conformal field theory. Since the dimer model considered in [28] imposes certain 
restrictions on the allowed dimer positions, it strictly speaking has not been proven 
to lie in the same universality class. 
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way for the matter to create short spin-boundaries is by squeezing individ- 
ual spatial slices to a minimal (cut-ofF) size, producing the long, stalk-like 
structure seen in Fig. 14. In the remainder of space-time, the spatial slices 
attain a macroscopic size. On this extended part of the universe, the scal- 
ing of the average spatial length in terms of the cosmological constant is 
anomalous, 

Wa~A^' (128) 

which was not the case for c < 1. Moreover, the critical Ising exponents, 
when measured on the extended region, are still given by the Onsager val- 
ues. We therefore have identified an analogue of the c = 1 barrier, at which 
a geometric phase transition takes place. However, its effect is much milder 
than in Euclidean gravity. Ising-life continues even beyond the barrier. This 
underlines yet again the drastic difference between 2d gravity with and 
without baby universes. 



6. Dynamically triangulated quantum gravity in d > 2 

6.1. GENERALIZATION TO HIGHER DIMENSIONS 

The general outline of the simplicial regularization of Euclidean quantum 
gravity presented above remains true in dimension d > 2, where the fun- 
damental building blocks are d-simplices with d+1 vertices. The partition 
function for dynamically triangulated pure gravity is given by 

Z = Y — 4— e-^^« , (129) 

where the sum is taken over a class of labelled d-dimensional triangulations 
with spherical topology. In the following, we will be mainly interested in 
the cases d = 3,4, where the discretized Einstein-Hilbert action (87) can 
be written equivalently as 

Seh = KdNd - kqNo, (130) 

using the Dehn-Sommerville relations. Since for simplicial manifolds and 
fixed volume A^^^;, the number of vertices has an upper limit Nq oc N^/d, the 
action (130) is bounded for finite volume. 

Note that in d > 2 it is not possible to realize a locally flat geometry 
by gluing together equilateral simplices. Only in 2d there is a regular tiling 
of flat space, in which each vertex is shared by six triangles. This happens 
because the dihedral angle 6d of each d-simplex contributing to the parallel 
transport around a (c?— 2)-simplex satisfies Oa = arccos 1/d, which for c? > 2 
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is not a rational fraction of 27r. However, this is a short-distance property 
which should become irrelevant in the scaling limit. 

In the simplest case of gravity without matter, there are two couplings: 
Kd, which is proportional to the cosmological constant, and kq, related to 
the gravitational constant. Both quantities are dimensionless in the lat- 
tice formulation, but will get their dimension back through "dimensional 
transmutation" in the scaling limit, if it exists. 

The gravitational state sum is now 

Z = Y,e-^''''ZviKo), (131) 

V 

where V = N^, and Zy{hq) is the partition function at fixed vohimc. In 
order for the sum (131) to be well defined, Zy{ko) must be exponentially 
bounded for large V, 

Zv(«o)~e«S(«o)V'^V'(/^o). (132) 

The existence of such a bound has been shown [24, 25], but the form of the 
subleading terms is at present unknown. A power behaviour of zy leads to a 
situation similar to (102), but there could in principle also be a logarithmic 
dependence logzy ~ 1/", with a < 1. For systems with a fixed volume 
one can replace Zy by zy when averaging over configurations. Note that 
higher-dimensional gravity with the action (130) resembles a spin system in 
2d, with the role of the matter coupling /3 played by kq- The cosmological 
constant will undergo a renormalization according to 

Z = ^e-('^<^-'^d)V(ACo), (133) 

V 

and — k^) gives an estimate of (V),^^- 

To understand the critical properties of the system defined by (133), 
one must investigate the limit V —>■ oo. For large positive kq; configurations 
with many vertices will be favoured. Comparing with the 2d case, we expect 
this limit to correspond to a branchcd-polymcr phase. Conversely, for large 
negative kq, the system will move towards lower iVo, leading to "crumpled" 
configurations. Note that unlike in 2d, Nq and V are independent, and Nq 
can be very small, for example, it is easy to find triangulations such that 
A^Q ~ V"" with a < 1. The fraction p{kq) = d{NQ)y /V can be used as an 
order parameter of the theory, and is related to the average local curvature. 
It may vary between and 1, corresponding to the two extreme geometric 
phases. In the limit p = (kq small), where the geometry is "crumpled", 
there are only a few vertices and one can move from one to another in only 
a few steps. The Hausdorff dimension of such geometries is large. In the 
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limit p = 1 (large kq), the number of vertices for a given V is maximal, 
leading to branched-polymer configurations with Hausdorff dimension two. 

Neither of these two limits looks like a promising candidate for a scaling 
d-dimensional quantum theory. Let us compare this situation with that of 
a 2d spin system. The two extreme phases correspond to the high- and low- 
temperature regions of the spin model, where also the spin system has no 
continuum limit. In the high-temperature limit, the spins fluctuate totally 
randomly, whereas in the low-temperature limit they hardly fluctuate at 
all. However, there is a critical temperature Tc where the spin fluctuations 
become critical and long-ranged. At Tc, one can take a continuum limit and 
recover a c=l/2 conformal field theory. In a similar way one might hope to 
find a critical gravitational coupling kq, where a transition between the two 
phases of geometries takes place, and where more "reasonable" geometries, 
exhibiting long-range fluctuations, can be observed. Such a point could 
serve as a non-perturbative fixed point where a theory of quantum gravity 
can be defined. 

Although the problem of solving higher-dimensional gravity with the 
simple action (130) is a combinatorial problem which superficially looks 
quite similar to the 2d case, wc still have no hints of how to find an exact 
solution. There is a mean-field analysis [25, 26] which gives a surprisingly 
precise description of the phase transition (see also [27]), but most of our 
understanding of Euclidean dynamical triangulations in d > 2 comes from 
performing Monte Carlo simulations, along the lines described in Section 
5. Also in these cases a set of ergodic moves is known [46]. Recall that 
in 2d, one possible set was given by the (p, g)-moves, for integer p,q > 
and such that p + q = d + 2. During a (p, g)-move, a subcomplex of p 
simplices is replaced by one of q simplices with the same (d— l)-dimensional 
boundary. This set can be generalized to higher dimensions [64, 65], using 
the observation that by gluing the two subcomplexes from before and after 
the move, one obtains the boundary of a ((i-)-l)-simplex. The boundary can 
be divided into two connected parts (a part X and its complement X) in a 
finite number of ways. Each move corresponds to replacing some X ^ X . 
In two dimensions, the boundary is that of a tetrahedron, made of four 
triangles. It can be divided either into two pairs of two adjacent triangles 
or into a single triangle and a set of three triangles sharing a vertex of order 
three. 

Performing an analogous construction in 3d, one obtains (1,4)- and 
(2, 3)-moves, together with their inverses (4,1) and (3,2). The first move 
inserts a new vertex of order four at the centre of a tetrahedron, thereby 
creating four new tetrahedra. A (2, 3)-move takes two tetrahedra which 
share a triangle, removes this triangle and replaces the two tetrahedra by 
three tetrahedra, whose common axis (a link of order three) is dual to the 
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Figure 15. The elementary ergodic moves in d = 3. There are four different moves since 
the moves are different from their inverses. 

removed triangle. Both types of moves are shown in Fig. 15. A similar 
construction in 4d leads to a set of five {p, g)-moves. 

These finite sets of moves are ergodic, that is, by repeating local moves, 
all configurations in the ensemble can be reached. Unlike in 2d, there are 
no subsets of moves which are ergodic in the set of triangulations of fixed 
volume. In 3d, all (p, g)-moves change the volume. If one wants to perform 
numerical simulations at some fixed V, a modification of the procedure used 
in the simulations becomes necessary. 

A simple way of forcing the volume to lie close to a given V is to use a 
modified action 

S'eh = i^dNd - KoNo + e\Nd - (134) 

in the Monte Carlo simulation, typically with (5 = 1 or 2, and a free param- 
eter e. The Monte Carlo process now takes place in the set of all triangu- 
lations, with fluctuating volume. If we succeed in tuning the cosmological 
constant to criticality, Kd ~ k'^{ko), the additional term in the action will 
force the volume to fluctuate around N^, = V. The amplitude of the fluc- 
tuations will depend on the free parameter e, which should be sufficiently 
small to permit big changes in the geometry. We can then simply collect all 
configurations with = V, and use them in measurements, provided they 
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are separated by a number of steps much bigger than the autocorrelation 
time. 

This method seems straightforward, but it leads to the following com- 
plication. Since our set of moves is ergodic only in the ensemble of simplicial 
manifolds of fixed topology and of arbitrary volume Nd, restricting the pos- 
sible volume fluctuations could make some parts of the configuration space 
unreachable. For manifolds which are not algorithmically recognizable it can 
be proven [66] (see also [8]) that the number of steps needed to connect two 
arbitrary triangulations of identical volume cannot be bounded by a 
computable fiuiction /(A^d) (a function which can itself be computed by a 
finite algorithm). This situation cannot arise when one has a finite set of 
moves which is ergodic for fixed volume (or even bounded by a computable 
function), since in this case one can just make a list of all triangulations, 
which will have a finite, computable length. We conclude that 

• ifa manifold is not algorithmically recognizable there will be "volume 
barriers" in the sense explained above, effectively separating parts of 
the configuration space. 

• In a dimension d where algorithmically unrecognizable manifolds exist, 
one cannot have a single, finite set of volume-preserving moves which 
at the same time is general, that is, which works for any manifold of 
the given dimension, and which is ergodic for all manifolds. 

In d = 2 all manifolds are algorithmically recognizable, and the flip move 
(introduced in Section 5.1) is a volume-preserving move which is ergodic for 
an arbitrary manifold. It is not known whether there exist three-dimensional 
manifolds which are not algorithmically recognizable. It has recently been 
proven that the three-sphere is algorithmically recognizable [67], but it is 
not known whether the same is true for 5^. For d > 3 there exist algorith- 
mically unrecognizable manifolds, and for dimension d > A not even the 
d-spherc is algorithmically recognizable [68]. 

In the 4d numerical simulations we use triangulations with S'^-topology, 
and the reliability of the procedure outlined above could depend on the 
manifold being algorithmically recognizable. If it is not, one might worry 
that there are configurations, or even large regions which are separated from 
the rest of the configuration space by very high "volume barriers". This 
would imply that in order to reach such a configuration by a successive 
application of the elementary moves, one would be forced to go through 
very large intermediate volumes. 

Attempts have been made to detect such barriers in numerical simula- 
tions on S^, but without success [69]. Unfortunately, this cannot be taken 
as an indication that is algorithmically recognizable because even for 
the five-sphere S^, which is known to not have this property, such barriers 
have not been seen [70]. There are two ways of interpreting the iS^-result. 
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Either the Monte Carlo simulation explores only a fraction of the config- 
uration space (in some sense it does anyway, since one always deals with 
a finite sample of "typical" configurations), or the configurations "beyond 
the barrier" which are hard to reach constitute only a very small fraction 
of the configuration space. For S"^, there need not be any problems in case 
it is algorithmically recognizable. If it is not, we hopefully are in a situ- 
ation where the "unreachable" configurations form only a small subset of 
configurations that becomes negligible in the infinite- volume limit. 

Just as in 2d, the local moves can be supplemented by minbu surgery 
moves, using a suitably generalized concept of (minimal) baby universes. 
They are again defined as parts of the triangulation which are connected 
to the parent universe by (d — l)-dimensional minimal necks, which topo- 
logically are boundaries of c?-simplices. The baby universes again play an 
important role, especially in the large-Ko phase, and using an improved 
algorithm has been essential in studying the phase structure of the four- 
dimensional theory. 

The observables described in 2d also exist for gravity in higher dimen- 
sions. The minbu distribution can be measured and used to determine ^str 
(in case the behaviour of zv{ko) is power-like), and the scaling properties 
of the two-point function Gv{r) yield information about the Hausdorff di- 
mension. Matter degrees of freedom can be included without problems, the 
simplest being spins and massless scalar fields. The coupling of gauge fields 
to 4d dynamical triangulations will be discussed in the next section. - In 
all these cases we can employ finite-size scaling and measure the matter 
correlators. 

6.2. NUMERICAL RESULTS IN HIGHER DIMENSIONS 

The algorithm outlined in the previous section was first applied to pure 
3d gravity [71]. As predicted, both a branched-polymer and a crumpled 
phase were discovered, but the transition between them turned out to be a 
very strong first-order phase transition. This manifests itself in a very large 
hysteresis in the order parameter p{kq), when changing kq across the phase 
transition point. The numerical effort spent on the 3d case was relatively 
small. Since the phase transition is first-order, it cannot give rise to a scaling 
limit. The phase structure remained unchanged also when spin fields were 
added [72]. 

Soon afterwards a numerical algorithm for 4d pure gravity was con- 
structed [73] and its phase structure studied extensively. It looked similar 
to the one in 3d, with two distinct phases. However, this time the phase 
transition was much softer, and for a long time was believed to be contin- 
uous. This exciting possibility led to many numerical explorations of the 
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properties of the two phases and the nature of the phase transition. It was 
shown that for large kq one obtains a generic branched-polymer phase with 
Istr = 1/2 and = 2. Minbu surgery moves were essential in generating 
independent configurations in this phase. The structure of the crumpled 
phase seemed rather unusual: the behaviour of zv{ko) was no longer dom- 
inated by a power of V, and the scaling of Gy(r) suggested that dh ^ oo. 
(This implies that the average linear extension of the system is almost in- 
dependent of V.) It was then understood that these pathologies are related 
to the appearance of exactly two "singular" vertices of very large order 
~ V. A finite fraction of the simplices of the triangulations shares these 
two vertices, giving rise to a "volume condensation". The two points are 
always joined by a link, whose order also grows with the volume [74]. 

The appearance of these structures is well described both by the mean- 
field approach mentioned earlier [26], and by a phcnomcnological, so-called 
balls-in-boxes model [27] . This model represents the gravitational system as 
a collection of A^o boxes, where the z'th box contains qi balls, analogous to 
the order o(i) of a vertex i in the siniplicial complex. Note that the vertex 
orders o{i) on a triangulation satisfy 



In the simplified model, the qi are required to satisfy the same constraint. 
The local probability of finding q balls in a box is given by p{q) ~ q^^ , 
where /3 is a free parameter. The model can be solved exactly in the ther- 
modynamic limit, and exhibits two phases: one in which the balls are dis- 
tributed randomly in all boxes, and a "crumpled" phase at large /3 where 
a finite fraction of the balls appears in one box. Note that this model has 
no concept of neighbouring boxes, and there are no correlations among the 
qi (other than the relation (135)). 

Although the balls-in-boxes model can hardly claim to reveal much 
about the geometric structure of quantum gravity, it nevertheless seems 
to capture some important features of the situation observed in the nu- 
merical simulations. In particular, one might wonder whether also in the 
full simplicial-gravity model there are no correlations between the vertex 
orders. This conjecture can be tested by investigating the corresponding 
gravitational two-point correlator for vertices separated by a (geodesic) 
distance r. Such a correlator is similar to a matter-matter correlator with 
non-vanishing expectation values for the matter fields. (We will not discuss 
how to define correctly its connected part.) Preliminary measurements in- 
dicate that the correlations do not extend beyond a few lattice spacings, 
that is. they are of the order of the cut-off [77] (see also [78] for an earlier 
report). Similarly short-ranged correlations have also been observed in the 
2d case [61]. 
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Figure 1 6. The double peak structure at the phase transion point as a function of the 
the number of vertices, TVo, for a fixed number of four-simplices, iV4 = 32000. 

There is by now overwhelming numerical evidence that the weak- and 
strong coupling phases are separated by a weak first-order transition [79]. 
The determination of the order turned out to be difficult, since no hysteresis 
was observed for finite systems. The best way of identifying the nature of 
the phase transition is by measuring the distribution of Nq in a sample of 
configurations exactly at the pseudo-critical point. The distribution has two 
peaks, whose separation grows with the volume. This is shown in Fig. 16. 
The effect becomes only apparent for sufficiently large systems (of 32000 
and 64000 simpliccs); for smaller volumes the separation cannot be seen. 

From the point of view of quantum gravity, this is of course a disap- 
pointing result, because there is no new continuum physics associated with 
the phase transition. The weakness of the transition could be an indication 
that a continuous phase transition can be reached by enlarging the space 
of coupling constants, that is, by extending the action Seh- Obvious can- 
didates would be terms containing higher powers of the curvature, which 
have previously been used in quantum Regge calculus [80]. 

In both 3d and 4d simplicial gravity, the phases observed on either side 
of the transition are pathological, in that they do not seem to correspond 
to an effective quantum geometry of a sensible dimensionality. As already 
mentioned we find either branched polymers (locally one-dimensional ob- 
jects) or a condensation of the volume around selected points (yielding an 
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object of arbitrarily large dimension). It is surprising that the rich geomet- 
ric structure observed in two dimensions has all but disappeared. There are 
potential cTircs for this, although there is no guarantee that they will lead 
to the desired quantum theory. 

Within the framework of dynamical triangulations, there are basically 
two ways of "changing the measure" , in order to improve this state of af- 
fairs. One is by using a modified action, but leaving the sum over geometries 
untouched. Gravity with a squared scalar-curvature term was investigated 
in [81], but the numerical results indicated that the phase structure re- 
mained unchanged. Note that also in d = 2, the addition of such terms 
does not change the universality class of pure gravity, as has been shown 
both numerically (see, for example, [82]) and analytically [83]. The other 
possibility is to change the sum to a different class of triangulations. A 
physically motivated example of this kind is the imposition of a causality 
constraint on Lorentzian geometries. It does lead to a change of univer- 
sality class in d = 2, as we have already seen. There are indications that 
also in d > 2 a causality condition acts as an effective regulator for the 
quantum geometry [84]. Numerical investigations of Lorentzian models in 
higher dimensions are under way. 

Another logical possibility is that only the matter-coupled theory pos- 
sesses a well-defined scaling limit. Recall that in two dimensions, there is a 
rather narrow window of physically interesting Euclidean continuum the- 
ories, namely, Liouville gravity interacting with conformal fields of charge 
< c < 1. We simply may have been lucky that pure gravity lies within 
this range. 

Matter interactions could in principle suppress the formation of baby 
universes in higher dimensions, and the inclusion of gauge (vector) fields 
was suggested as a promising scenario for testing this conjecture [85]. The 
simplest model of this kind is the non-compact 4d Abelian gauge model 
formulated in [75]. Like on a regular lattice, the gauge fields Aij = —Aji 
are associated with oriented lattice links A discretized version P of 

the field strength tensor Ff^^i, can be expressed in terms of "plaquettes" 
around (oriented) triangles {ijk}, 



(136) 



The action was taken to be 



S, 



'gauge — 



(137) 



{ijk} 



with o{{ijk}) denoting the order of the triangle {ijk}. This factor was in- 
cluded to guarantee that whenever is constant for all triangles, Sgauge 
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is proportional to the volume of the manifold. The action has a sim- 
ple gauge symmetry and is relatively easy to use in numerical simulations, 
since it is Gaussian in the vector fields. The number Ug of gauge fields 
can be varied, and for Ug = 3 a new "crinkled" phase with 7^^^ < was 
observed, replacing the branched-polymer phase. Moreover, the phase tran- 
sition seemed to be continuous rather than first-order. 

However, it was soon realized that this new phase can be obtained in a 
model without gauge fields, by adding a geometric term 

Sa = ^ E ^ogo{{ijk}) (138) 

{ijk} 

to the pure-gravity action [76]. Since Sa is an ultra- local measure term, 
one would not really expect it to lead to new continuum physics. 

The origin of this term can be traced to a somewhat unfortunate choice 
of gauge action. At first sight, it seems natural to associate the gauge fields 
with the links of the triangulation, analogous to what is done on regular, 
flat lattices. This implies that any charged fields which couple to the gauge 
fields will be located at the vertices. However, note that scalar fields coupled 
to gravity in all investigations so far have been placed at the centres of the 
four-simplexes, to ensure that the number of field degrees of freedom grows 
proportional to the volume. In the case of fluctuating geometry, this need 
not be the case if the scalar fields are located at the vertices, since Nq may 
grow much slower than A^4. 

Thus, if our primary aim is to compare the effects of gauge fields on 
geometry with those of the previously studied scalar fields, the complex- 
valued charged fields should be associated with the four-simplexes of the 
triangulation, and the gauge fields with the links of the dual lattice. The 
pure-gauge action is then a sum over the dual plaquettes, which are the 
natural duals of the triangles in the original triangulation (see [77] for 
details) . 

The two models with the gauge potentials A on the links and on the 
dual links are dual to each other. The theory defined by (137) (on the 
triangulation) is equivalent to the gauge theory on the dual lattice, if we 
add by hand a term (138) to the natural gauge action on the dual lattice 
[77]. Monte Carlo simulations with gauge fields on the dual links without 
the term (138), as advocated above, have not detected any effect of the 
gauge fields on the geometry, in line with previous results for scalar fields 
and Ising spins coupled to gravity. 

Putting a more positive spin on this result, the inclusion of the term 
(138) is the only example of a modification of the action which has led 
to a genuine change of the geometry. Although it has no obvious physical 
interpretation and although the crinkled phase is probably not related to 
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interesting continuum physics, it does give us some hope of finding a physi- 
cally well-motivated term which does change the geometry in an interesting 
way, allowing us to define a continuum theory of gravity. 

7. Outlook 

In these notes wc have described in some detail a non-perturbative path- 
integral approach to quantum gravity, with and without matter. This "dy- 
namical triangulations" method uses an intermediate regularization, in 
which the space of all space-time geometries is approximated by a set of 
simplicial complexes with certain edge length assignments, in the spirit of 
Regge calculus. For a finite volume (that is, for a complex with a finite 
number of simplices) and fixed space-time topology, this method yields 
well-defined, convergent partition functions and propagators. The aim in 
this approach is to understand the critical behaviour of these statistical 
systems of fluctuating geometry, and investigate their physical properties 
in the continuum limit, if it exists. 

We believe that this ansatz has a number of virtues that make further 
research worthwhile. Firstly, it works perfectly in two dimensions, leading 
to a host of analytical and numerical results. There, it agrees with other 
discrete or continuum formulations, whenever they can be compared, and in 
some cases turns out to be even more powerful. This gives us considerable 
confidence in the validity of the method. Moreover, the two-dimensional 
construction has served as a useful laboratory for defining diffeomorphism- 
invariant versions of two-point functions and other field-theoretic observ- 
ables. Dynamical triangulations as a method can be applied in any dimen- 
sion d, but obviously will lead to very different theories, depending on d. 
Only in d = 4 do we expect to find local, interacting quantum-gravitational 
degrees of freedom. 

Secondly, as has been shown recently, there is also a Lorentzian version 
of the formulation, where additional causality restrictions are imposed on 
the space-time geometries, and a well-defined Wick rotation exists. The 
absence of a natural notion of "analytic continuation" from Euclidean to 
Lorentzian signature has for a long time stalled progress in discretized path- 
integral approaches to gravity. Surprisingly, in two space-time dimensions 
the Euclidean and Lorentzian continuum quantum gravity theories are very 
different. As a consequence of the causality constraint, the "baby- universe" 
structure of the Euclidean theory is absent in the Lorentzian sector. This 
leads to a "smoother" quantum geometry, and well-behaved matter cou- 
pling even beyond the c = 1 barrier. Such a behaviour is very desirable in 
dimension d > 2, where so far the dominance of "extreme" geometric phases 
in Euclidean dynamically triangulated models seems to have prevented the 
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existence of an interesting continuum limit. 

In higher dimensions, much work remains to be done to estabhsh whe- 
ther our method leads to a well-defined continuum quantum theory. In the 
Euclidean sector, a suitable change of the discretized measure or the inclu- 
sion of different types of matter may still lead to long-range correlations. 
In the Lorentzian sector, there is a genuine chance that the causality re- 
strictions will lead to a qualitatively different phase structure, including a 
second-order phase transition. The investigation of the Lorentzian models 
for d > 2 has only just begun. However, we do not have to start from 
scratch when trying to solve the theory, but can use the considerable ex- 
pertise gathered in the ETiclidean simulations. Obviously, quantum gravity 
will not be solved by numerical methods alone, but we have learned from 
the example of dynamical triangulations in 2d that a mathematically exact 
analysis and numerical simulations can work hand in hand. 

In addition, one can hope for analytical progress in solving the dis- 
cretized models (which in the case of quantum Reggc calculus has always 
been hampered by the presence of inequalities on the edge lengths, even 
in d = 2). In terms of dynamical triangulations, quantum gravity has been 
turned into a well-defined combinatorial problem, which has already been 
solved exactly in two dimensions. It is also likely that - analogous to the 
Euclidean case - mean-field methods can help to elucidate the Lorentzian 
phase structure in higher dimensions. This may have interesting ramifica- 
tions for the general study of spaces of pseudo-Riemannian geometries. 

In summary, the dynamical-triangulations approach is a promising tool 
for both analytical and numerical investigations of quantum gravity. It al- 
lows for a natural inclusion of matter, and is explicitly reparametrization- 
invariant. The path integral can be defined for both Euclidean and Lorentz- 
ian geometries. In the Lorentzian case, it involves summing over a class of 
Wick-rotatable space-times with a natural "time" -slicing, reminiscent of 
the geometric structures appearing in Hamiltonian quantization schemes. 
We hope that our construction of a Lorentzian path integral will serve as a 
bridge between the canonical and covariant formulations, benefitting both 
and improving our understanding of the nature of quantum gravity. 
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